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ON THE DEGREE OF APPROXIMATION TO A 
FUNCTION BY THE CESARO MEANS OF ITS 
FOURIER SERIES 
By T. M. FLETT (Liverpool) 


[Received 2 August 1955] 


1. Let f(x) be integrable in (—z, 7) and be periodic with period 27, and 
let x cs 
f(x) ~ 4ao+ (a, cosna+b, sinnx) = fa9+ ¥ p,(2). (1.1) 
1 1 


It is a classical result of 8S. Bernstein¢ that, if f(x) is of class Lip «, where 
0 <a« < 1, and o,(x) is the nth Fejér mean of the series (1.1), then 
o,,(x)—f(x) = O(n-*), (1.2) 
uniformly in x. This result continues to hold for the (C,*) means for 
any positive k. For the partial sums s,, however, we have only 
8,(x)—f(x) = O(n-“log n), (1.3) 


uniformly inx.t Salem and Zygmund (13) have shown that the logarithm 
in the relation (1.3) cannot be suppressed even if, in addition to the 
hypothesis that f belongs to Lipa, we suppose also that f is of bounded 
variation. They show, however, that, if f belongs to Lipa and is of 
monotonic type [i.e. there exists a constant C such that f(x)+ Cz is 
either non-decreasing or non-increasing in (—00,00)], then we do have 
8,(%)—f(x) = O(n-*), (1.4) 
uniformly in 2. 

There is an essential difference between the results (1.2) and (1.4) [or 
(1.3)] which becomes apparent only when we consider these relations for 
a particular value of x. For, while the property (1.4) [or (1.3)] depends 
on the behaviour of f in the entire interval (—7, 7), the property (1.2) 
depends only on the behaviour of fin the neighbourhood of the particular 
point x concerned. To prove the first statement it is enough to remark 
that, on the one hand, the property (1.4) implies that p,(x) = O(n-*),§ 
while on the other hand, if the only non-local hypothesis is the integra- 
bility of f, then p,(x) = o(1) and need not be of smaller order [for 

+ See e.g. Zygmund (14), 62, Ex. 7. 

t This result is due to Lebesgue. See e.g. Zygmund (14), 61, Ex. 6. 

§ This follows by a trivial modification of the well-known ‘limitation theorem’ 
for Cesaro summability [see e.g. Hardy (5), Theorem 46]. 
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82 T. M. FLETT 


¥ p,,(x)cos nt is the Fourier series of the function ${f(x+t)+f(x—t)—ay}]. 
The second statement is the case k = 1 of the following known theorem, 
in which o*%(x) denotes the nth (C,*) mean of the series (1.1), and 


P(t) = ¢,(t) = f(x+t)+f(x—t)—2f(z). 

THeoreM A. Let 0<a<10<i<a,k >a. If xis a point such 
thatt $,(t)| < At® (1.5) 
when 0 <t <5, then 

ok(x)—f(x) = O(n-*). (1.6) 

The result of Theorem A holds under weaker conditions than (1.5), the 
most general result of this type being that of Obrechkoff (12).t Wereturn 
to the question of these more general conditions at the end of this paper. 

The relation (1.6) is still a local property when k = a, but in this case 
the condition (1.5) is no longer sufficient for the truth of (1.6). We prove 
here the following result.§ 


THEOREM 1. Let 0 <a <1,0<8 <a. If xisa point such that 
t 
| |\dd,(u)| < At (1.7) 


0 


when 0 <.t <8, then 


ox(x)—f(x) = O(n-). 


The condition (1.7) of Theorem 1 cannot be replaced by (1.5), nor can 
we increase the order of the expression on the right of (1.7). More pre- 
cisely, we prove 


THEOREM 2. Let w(t) be a positive function, increasing to infinity with t 


as slowly as we please. Then for each positive a less than 1 there exists a 
function f, of bounded variation in (—7, 7), such that 


(i) $o(t) = Olt); 
(ii) [ \dog(u)| = Ofteu(t-)}; 
0 


(iii) an(0)—f(0) ~ O(n-*). 


+ Luse A(6, c,...) to denote a positive constant depending only on 6, ¢,...; A by 
itself will denote a positive absolute constant. These are not necessarily the same 
on any two occurrences. The constant implied by an O will in general depend 
on the various parameters concerned. 

t I am indebted to the referee for this reference. 

§ I am indebted to Professor J. E. Littlewood for the suggestion that the 
condition (1.7) might be relevant. 
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For k < « the relation (1.6) is no longer a local property, for it implies 
thatt p,(x) = O(n*-*), (1.8) 
and this is again smaller than need be the case if the only non-local 
hypothesis is the integrability of f. This is in fact the only reason why 
the validity of the relation (1.6) ceases to be a local property when k < «. 
For, if we take (1.8) as a hypothesis, we find that the relation (1.6) depends 
only on the behaviour of f in the neighbourhood of the particular z con- 
cerned. Corresponding to Theorem 1 we have 
THEOREM 3. Let O<a<1O<B<10<8i<2,k >a—f. Ifzx 
is a point such that 


p, (x) = O(n-*), 


t 
and dd (u)| < At* 
6 
when 0 t < 8, then 
ok(x)—f(x) = O(n-*). (1.9) 


The problem of localizing the property (1.9) is very similar to the 
problem of localizing summability (C,k) for negative k discussed by 
Bosanquet and Offord (4), the O(n-*) of (1.9) there being replaced by 
o(1). The fact that here we require only an O-relation, however, enables 
us to use somewhat simpler arguments. 

It is easy to see that, if f(2) belongs to Lip « and is of monotonic type, 
then the conditions of Theorem 3 are satisfied at every point x for every 
8 < 1. We therefore have the following corollary, which contains as the 
special case k = 0 the result of Salem and Zygmund mentioned earlier. 

THEOREM 3, CoroLiaRy. If f(x) belongs to Lipa, where 0 <a < 1, 
and is of monotonic type, and if k > a—1, then (1.9) holds, uniformly in x. 

There is a corresponding result for the case « = 1, and here we can 
include the case 8 l. 

TuHeoreM 4. Let 0 <8 <1,0<i<2,k>1-8. If xis a point 
such that p, (zt) = O(n-*), 

: 
and | dd,(u)| < At (1.10) 
} 0 
when 0 <t <3, then 
ok(x)—f(x) = O(n log n). 


The condition (1.10) is implied by the condition that f is of class Lip 1 


+ See.the third footnote on p. 81. 
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in the closed interval (2—8,2+-5).t We therefore have the following 
corollary. 

THEOREM 4, COROLLARY. Suppose that f is [integrable in (—7, 7) and] 
of class Lip 1 in the closed interval (a,b), and that a,,, b,, = O(n-*), where 
0<B <1. Then, fork > 1-8, 

ok(x)—f(x) = O(n log n), 
uniformly in any closed interval interior to (a,b). 

There is a corresponding result for the other Lipschitz classes, namely 

THEOREM 5. Suppose that f is integrable in (—7, 7) and of class Lip « 
in the closed interval (a,b), where 0 < « < 1, and that a,b, = O(n-), 
If0 <B <aandk > «—f, then 

ok(x)—f(x) = O(n-*), 
while, if B = «, s,(x)—f(x) = O(n-“log n), 
both results holding uniformly in any closed interval interior to (a, b).t 


We prove Theorems 1-5 in $§ 2-4. In § 5 we extend these results to 
the conjugate series, and in §§ 6-7 we consider generalizations of these 
theorems under weaker conditions. 

The substance of the remarks in § 6 was suggested by the referee, 


to whom I am also indebted for many helpful comments. 


2. Proof of Theorem 1. We have 


ok(x)—f(x) = : [ A(t) K*(t) dt, 
where 4(t) = ¢,(t) = f(zw+t)+f(~a—t)—2f(x), and K*(t) is the nth (C, k) 


mean of the series $+cost+-cos 2t+-.... Let (n;t) stand for 


{(n+4k+4)t— kr}, 


T'(in+k+1) ve 
P(k+1)P(n+1) 





and let B= 


+ That is to say, that |f(u+t)—f(u)| < A|t} whenever u+?t and w lie in the 
closed interval (a#—8, 2+8). 

t It has been pointed out to me by the referee that Theorem 4, Corollary, and 
Theorem 5 belong to a class of localization theorems, considered by Zygmund, 
in which conditions are imposed on a, and b, and the conclusion holds for an 
interval [see e.g. Zygmund (14), 305, Ex. 12]. There is thus a distinction between 
these results and Theorem 3, in which the condition bears only on p,(2) and the 
conclusion holds only for a single point. 
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It is known? that, for 0 <t < 2, 


K*(t)| < A(k)n, (2.1) 
while for n-! < t < mw we have 
Kk(t) = R(t)+ S(d), (2.2) 


a i) 
~ BE(2 sin $t)F+!’ 
and S(t)| < A(k)n-'t-?. (2.4) 


We prove now two lemmas from which the desired result follows. 


where R(t) 


Lemma 1. Jf d(t) is integrable, and if0<a<1,0<8<2,k >a, 
then { d(t)K*(t) dt = O(n-2). 
Lemma 2. 1f0 <a <10<i<2,k >a-—l, andif 
t 
| |\dd(u)| < At* 
0 
when 0 t <3, then 
{ d(t)K¥(t) dt = O(n-2). 


_— 
tho 
or 

— 


Lemma | is an tmmediate consequence of (2.2), for in (5,7) we have 
K*(t) = O(n-*). To prove Lemma 2 we writet 
3 lin ry 5 


( dK dt = ( dK*dt+ (dRdt+ [ dSdt=1,414+h. (2.6) 
0 0 lin In 
Since 4(0) 0, we have 
} t 
d(t) d(t)—(0)| = | | db| < | |dd| < Af?, 
‘9 0 


whence, by (2.1) and (2.4), 


l/n 
I, O(n | rai = O(n-*), 
0 
8 
fn | ees a = O(n-*). 


1)n 


on 
wo 


Next, let A(t) = [ R(w) du. 


+ See e.g. Hardy (5), 360. 
t The method used here is due to Hardy and Littlewood (6). 
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By the second mean-value theorem, 
A(t)| < A(k)n-*-1t-*-1, 


Hence 
3 5 3 3 
I,= | ¢Rdt=— [ gA’dt=— | pdA= —[GA},+ | Add. 
Ij/n ljn ln Ij/n 
The first term on the right is evidently O(n-*). Further, if ¢*(¢) is the 
total variation of ¢ in the interval (0,¢), then |¢*(t)| < At*, whence 
| 3 a. : 
| | Add| < | |A| dg* < A(k)n-* | t-*-1 dg* 


| 
| 
l/n lj/n I/n 


C) 
= A(k)n-* ffeeagep +041) | t-k-24* a 
1jn 
: 
= O(n *)+0(n- [ t2-*-2 a 
1/n 
= Ofa-*), 
as required. 
Theorem 1 now follows from Lemmas | and 2. From Lemma 2 we 
obtain also the following result, which, although contained in Theorem 3, 
gives a simpler proof of the corollary of that theorem. 


THEOREM 6. Let f(x) be of bounded variation in (—7, 7), and periodic 
with period 27, and let0 <a <1,0<8<2,k >a—l. If xisa point 
such that t 

[ \dd,,(w)| < At 
0 
when 0 <t <8, then 
ok(x)—f(x) = O(n-*). 

To deduce this from Lemma 2, we have only to observe that, if 
(2.7) holds for some 5 < 7, and if f is of bounded variation in (—z, 7), 
then (2.7) holds with 6 = z. 


3. Proof of Theorems 3, 4, 5. We require the following lemma. 
Lemma 3.¢ Let —1 < k < 1, let 
tlh a — 4¢)*- 
where 18 = (2sin 38)-*-!, and suppose that 
pa(x) = O(n-#), (3.1) 


+ This is similar to a result of Bosanquet and Offord (4), but the proof is 
somewhat simpler. 
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where 0 <B <1. Then 


( d,(t)x(t)sin(n; t) dt = O(n -B), 


0 
We have $,(t) +2 x)cos vt, 


where Po = 4—2f(z). 
Hence [since y(¢)sin(n;t) is of bounded variation in (0, )}] the integral 
on the left of (3.2) is equal to 
Po ( x(t)sin(n; t) dt +-2> p, ( x(t)sin(n; t)cos vt dt = &, 
0 a 
Writing p_, = p, we get immediately 
==> p, [ x(fsin(n—v; 8) de. 
2°" 3 

Now, if n—v+43k+}4 + 0, we have 

r , , -0O8 ——p? f id —_ 

y(t)sin(n —v;t) dt : x(t)cos(n —v; t) ™ P ‘(t)cos(n—v; t) dt. 
‘ n—v+$k+} ], n—v+tk+4 
0 0 (3.4) 


a 








Since x(0) cos(n—v;7) = 0, the integrated term vanishes. The second 


term on the right of (3.4) is equal to 


; os }t 
or Fi cos(? v;t) dt T (k . 1) Sain deers CO” f) a 
0 5 
and, by the second mean-value theorem, this does not exceed in absolute 
value 


Hence 
—1 : x 
> py | x(é)sin(n—v; t) dt = 0} > v-A(n+v)-*| 
T 


0 


n+1 


= O(n- ‘Zs 8) 4 O( ¥ »-2-*) 
‘B), 


= O(n 


+ See e.g. Zygmund (14), 91 
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Similarly 
n—l n—1 (in) n—1 ’ 
>= 0} * v-8(n—v)-*| O\n-* > yBl +O\n-* Y (n—v)-* 
. : 1 (in}+1 
= O(n-!-8)+ O(n-8) = O(n-4), 
and 


al <= .-8.. -2| _ pln-8 < .-2) — sal 
a O12) B(v—n—1)-*} = O(n p> 2) = O(n-*). 
Since the remaining terms on the right of (3.3) are trivially O(n-*), this 
proves the lemma. 


Consider now the proof of Theorem 3. By Theorem A, we may suppose 
thata—8 <k <a. By Lemma 3, 
. : 4 
nfok(x)—f(x)} = | KE dt + | dRdt + | dSdt 
i) 5 8 
: 
Bk | th sin(n; t) dt +O(n-*-8). (3.5) 


0 


; : 
-” [ $Kk dt + | gS dt — 
3 5 


By (2.4), the second integral on the right of (3.5) is O(n-!) = O(n-*), 
while the first integral is O(n-*), by Lemma 2. It remains only to show 
that 


5 
i | , 
EE th sin(n:t) dt = O(n-*). 


0 

By hypothesis, ¢ is of bounded variation in (0,5), whence ¢ = ¢’—¢’, 
where ¢’ and ¢” are non-negative and non-decreasing in (0,5). By the 
second mean-value theorem, 

) 

{ 1d’ (t)sin(n; t) dt = O(n-), 

0 
and similarly for 4”. Hence 

; 

a th sin(n; t) dt = O(n-*-!) = O(n-*), 


ne 
0 


as required. 


To prove Theorem 4, we have only to make the obvious changes in 
the proof of Theorem 3 above, the expressions on the right of (2.5) and 
(3.2) now being replaced by O(n-' log n). 











APPROXIMATION BY CESARO MEANS 89 


To prove Theorem 5 we require 


Lemma 4. Suppose that f is integrable in (—7, 7), and of class Lip « in 
the closed interval (x—8,x+-8), where 0 <a <1. Then 
3 
Se , 
(1) EE th _(t)sin(n;t) dt = O(n-*) (k > 0); 


_ 


: - bl 
(ii) | b,(t)KR(t) dt = fO(n-") (k > 0), 


~— (O(n “*logn) (k= 0). 

The proof of both results follows almost exactly Lebesgue’s proof of 
(1.3),+ and we therefore omit it. Theorem 5 is now an immediate conse- 
quence of Lemma 4 and (3.5). 


4. Proof of Theorem 2.t For each positive integer v let n, = 2*, let 


+ 


m, = n,+4a+-4, and let p, be the smallest integer exceeding w(2m,_,/a7). 


y 


We may evidently suppose w increasing so slowly that 
a, < are (4.1) 


for all vy. Let J, , denote the interval 


(da +p—l1)z/m, <t< (4a+/p)z m,. 
For all large enough v, say for v > vg, we have p, > 2 and 
(a+ p,)7/m, < fan/m,_, 


> 


[the second inequality being a consequence of (4.1)], so that for v > vp 


the intervals J, ,, where p = 1, 2,..., p,, do not overlap the intervals J,_, .. 
Finally, let f be defined in (0,7) as the function which is equal to 
(—1)?-1(2 sin }#)* in J, , (v = v9, v9 +1...., p = 1, 2,..., p,), and to 0 else- 


where, and let f(—t)= (f(t). Evidently |f(t)} < Ajt}* for all ¢ in 
(—7W, W). 
The first part of the argument of Lemma 2 shows now that 
7 


; 9 2 sin(n; t) - 
o2(0)—f(0) = =f f(t)K2(t) dt = om: | IO Sanya @ +O(n-*), 
0 7 


+ Zygmund (14), loc. cit. 
t The function constructed here is similar to one used by Izumi (8, 9), and 
Izumi and Sunouchi (10), for various purposes in the theory of Fourier series. 
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where (n;t) stands for {(n-+4a+}$)t—4a7m} and yn = fam/(n+4a+4). 
is therefore enough to prove 


(ii) the total variation of f in (0,t) is Ofts(1/t)}; 


iho we sin(n; ft) 
(ill) t)_—__——_——. dt + O(1). 
4 I @sin 1) . ”) 
7 
Consider first (iii)’. We take n = n,, where » > v9, and for any 
interval J = (a,b) we write 
; a 
(2,9) (2 sin 1 
In the intervals J, , for p = 1, 2,..., p,, the function f has the sign of 
sin(n,,:t), whence 
* sin(n,,; t) * |sin(n,,:t) 


H—— > en OE 
Pain i : t 


J (I 


eal 


. 9 


a | sin(n;t)| dt = (Ja+p)= ai 


Bans Tine 
Since f = 0 when 

+p, rm, <t <S 3a 

we have Pu 

J (Saz/m,,, }a77/m,_4) +A 2 (etl) > A log p,,. 


On the other hand, for any interval J, , we have 


P om, st) - Am, 


J (I, ,) 2 ..- <n 
np 2 sin Wt m,(4a+p—1) 


Lip 
by the second mean-value theorem, whence, if v < p, 


py 
J (ham/m,, tam/m,_,) < Amz*m, > (4a+p—1)"! 
p=1 


< A(a)my*m, log p, < A(x)n7*n, log n,,. 
Thus 
p-l 


J(ham/m,-1,7) < A(a)nztlogn, > 2” 
v=Vo 


< A(a)nz* log >? 
< A(a)n;*logn, 2°* = A(a)nz*logn, = 0(1) 


Til 
(4.3) 
as 100. Combining (4.2) and (4.3) we obtain (iii)’. 
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[t remains only to prove (ii). The total variation of f in the interval 
(0, t)is obviously of greatest order when tis one of the points ($a+-p,)7/m,. 
Now the total variation of f in either of the intervals 
san m <t< sam m,_4 


han/m, <t < (40+ p,)7/m,, 


Py 
is O(m>* ¥ p*) = O(n>* p}**). 
I 


If now t = (}a+-p,,)7/m,, where p > vo, the total variation of f in (0, t) is 


r 


Ol & n>* p}**) = O(nz* Pu’ 3 
v= pe 
for, by (4.1), 
> => 0% < D> 2-@ < Alan". 
v=pt+l v=p v=2H 


Since ¢ is of exact order n> the total variation of f in (0, ¢) is 


pe Pw 
O(t*p,,) = Oft*w(2m,,_ x7)} = Oft%w(1/t)}. 
This proves (ii), and completes the proof of Theorem 2. 
5. Results similar to Theorems 1-6 hold for the conjugate series of f, 


> (—b,, cosnx+a, sinnx) = > p,(z), (5.1) 
i I 


but here there is so difference between the cases a < 1 and a = 1. 
Let &(x) be the nth (C,k) mean of the series (5.1), let 
u(u) = b(u) = f(z+u)—f(x—u), 
and let f be the conjugate arf i.e. 
7 


f(x) = — ys,(u)cot fu du.t 
Then we have, for instance, 


THEOREM 7. Let 0 - C1O<Bp<10<i<2a,k>a-f. Ifz 


is a point such that 


p(x) = O(n-*), 

t 

and | |dyb,(u)| < At 
0 

when 0 <t <8, then 


zk (x)—f (x) = O(n). (5.3) 


+ Our conditions imply in particular that the integral on the right of (5.2) is 
absolutely convergent. 





92 T. M. FLETT 

If f is of class Lip 1, the conditions of Theorem 7 are satisfied at every 
point x for every 8 < 1, and we therefore have the following corollary 
[ Alexits (1), Zygmund (15)}. 

Corotiary.t If fis of class Lip 1, then (5.3) holds for k > 0, uniformly 
in x. 

The only novelty here is the inclusion of the case « = 1, and for this 
we require an estimate for the corresponding kernel, which does not seem 
to have appeared in print before. Thus 

7 
' 2 i 
F(x) = —— | w(t) L(t) dt, 


7 


0 


where L*(t) is the nth (C, &) mean of the series 0+sin¢+-sin 2t+-.... For 


vai qe, L¥(t)| < A(k)n, (5.4) 
while for n-! < t < x [indeed, for 0 << t < 2, ifk < 2] 
L¥(t) = }cot #+U(t)+Vit), 
cos(n; t) 
E*(2 sin $t)*+)’ 
andt V(t)| < A(k)n~*t-. (5.5) 

The inequality (5.4) is known. The inequality (5.5) may be deduced 
from the corresponding result for A*(t), or may be proved directly by 
the method used by Hardy (5, 360) and Kogbetliantz (11) to estimate 
K*(t). 

Using these inequalities in place of those of (2.1)-(2.4), we have only 
to make the obvious changes in Lemmas 1-3. It should be noted that, 
in virtue of (5.5), the integral J, in (2.6), which produces a logarithm in 
the Fourier series case when « = 1, is now innocuous. 


where U() = 


6. We consider next a more general condition for the validity of 


Theorem 3. For t > 0 write 
§ 
®,(t) = 2 (t—uy—'d(u)du (r > 0), 
I(r), 


0 
®,(t) = d(¢), ®,(t) = O1.,(t) (-l<r< 0). 


+ The result of the corollary may also be obtained more simply from the 
analogue for the conjugate series of Theorem 6. 

t Estimates with A(k)n™t? on the right of (5.5) are known. See e.g. Bhatnagar 
(2), Zygmund (14), 49. 
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It is known {Obrechkoff (12)] that, if for some r > 0 and 0 < a < 1 


t 


} u-"|D(u)|du << AP** (OS t <b <7), (6.1) 
0 

then ok(x)—f(x) = O(n-*) 

for k > ar? 


The gap between the condition (6.1) [which is a generalization of the 
condition (1.5) of Theorem A] and the condition 


( idd(u)|.< At* (0<t<8<72) (6.2) 


0 
of Theorems 1 and 3 can be bridged by a condition of the form 


t 
[ u-r'd®, .(u)| < AP? (OS t<b <2). (6.3) 
r 
A condition of this type [with a = 0] was first introduced by Bosanquet 
(3). The arguments of that paper show that, if (6.2) holds, then (6.3) holds 
and ®,,,(+0) = 0, for every r > — 1.{ The condition obtained by 
putting 7 | im (6.3) is actually weaker than (6.2), but implies (6.2) if 
in addition 4(t) = O(t*). Bosanquet’s arguments show also that, if 
(6.3) holds for a particular r = ry, > —1, and if ®,,,(+0) = 0, then 
(6.1) holds for every r > ry [so that, in the theorem of Obrechkoff’s 
quoted above, (6.1) can be replaced by (6.3) together with the condition 
®, (+0) 0}. 
The result corresponding to Theorem 3 with (6.2) replaced by (6.3) is 


as follows: 


THEOREM 8. Let O<a<10<B<1, —-l<r<9Q, 
> a—B, k > a+r, and let x be a point such that 


(i) p, (x) = O(n-*); 


(ii) ®,.(+0) = 0 


t 
+ The integral here is to be interpreted in the first instance as lim { , where 


e—0 € 
it is assumed that ®,., is of bounded variation in every interval (e,¢). When 
(6.3) holds, however, it can be shown that the integral on the left of (6.3) exists 
as an ordinary Lebesgue—Stieltjes integral [see Bosanquet (3), 114, footnote]. 
+ We use here the fact that (0) = 0. 


+ 
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and ‘ 


{ u-rid®, ..(u)| < AP+® (0 <t <8); 


0 


t 


(iii) { d(u) du = O(t'+), 


0 
Then ok(x)—f(x) == O(n-*). 
With the notation of § 2, we have 
ln 7 7 
n{ok(x)—f(x)} = [ Kk dt + [ dRadt + [ oS dt. 
0 ln ln 
By means of the condition (iii) we can show easily that the first and last 
integrals on the right are of the required order.¢ By (i) and Lemma 3, 
we have further 
. ; ; 
dR dt — oR dt — | td(t)sin(n; t) dt + O(n-*-8) 
ln lin 0 
. - 
— dR at — td(t)sin(n; t) dt +O(n-*-8). (6.5) 


me 
ln l/n 


The proof can now be completed by applying to each of the integrals 
on the right of (6.5) an argument used by Bosanquet (3, Theorem 1), 
using the condition (ii). Since the details are similar to those in Bosan- 
quet’s paper, I omit them, 

No doubt a result similar to Theorem 8 holds for the conjugate series, 
but I have not pursued this point. 


7. We note in conclusion that, when k < a, the condition (iii) of 
Theorem 8 is necessary for the truth of (6.4). We have in fact 


THEOREMS Jf0<k <a < 1, and 
ok(x)—f(x) = O(n-*), 


t 
then | d(u) du = O(t'**). 
0 
This result is a companion to a theorem of Hyslop (7), which states 
that, if (7.1) holds fork > a > 0, where 0 < a < 1, then ®,(t) = O(t*") 
for r > k—a+1. It is also a companion to a theorem of Hardy and 
Littlewood (6) which states that (7.2), with the O(f+*) on the right 


+ See Bosanquet and Offord (4) for a similar argument. 
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replaced by o(), is a consequence of the Cesaro summability of the 
series to any negative order. 


Theorem 9 is an immediate consequence of the following theorem. 


THEOREM 10. Let ok be the n-th Cesdro mean of order k of the series 


< 
» ‘ - 1 ‘ 
2 a,. If0<k 1, and for some 8 


ok —s = O(n-*), 
then, as t > 0-, 
Sa, 2™ _2 = oe). 
= nt 


The proof of Theorem 10 is similar to that of the result of Hardy and 
Littlewood quoted above, and we omit it. 
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ON THE LOWER ORDER OF ENTIRE 
FUNCTIONS DEFINED BY DIRICHLET 
SERIES 


By Q. I. RAHMAN (Aligarh) 
[Received 11 August 1955) 


1. ConstDER the Dirichlet series 


n 
no 


f(s) = Sa ,e* (Ayu >A,; A, > 0; limA, = w; s = o+it). 
1 


Let o, and oc, denote the abscissa of convergence and the abscissa of 
absolute convergence, respectively, of f(s). If 0, = 0, f(s) defines an 
entire function. I shall suppose throughout that o, = o, = ©. 


Let u(c) = maxia,, |e», M(c) = max |f(o+it)). 
n2>1 —ao<t<a 


It is known [(1) 68] that u(o) < M(o). 


On the other hand, if 


; log 
lim sup—-— = D < @, 


‘ 


nn n 


then [(1) 68] 
M(c) < p(ot+D+e) (€ > 90; o > ofe)). 


If (2) holds, then [(1) 69-71] 
log aA loglog M(c) 


; A , (aes 
lim sup “—_-_* = p = limsup ————_ (0<p<o). (4) 
n> og ay ox Co 


Under the same condition, a similar result for the lower order A, namely 


loglog M(c) (5) 
Co 


lim inf A, logan = A = liminf 


nx og a n| -1 o-oo 


does not always hold. In fact, for 


exp(e**)+exp(e*) = arerrem( a —~ zi) 


log A 


or — ] 
? 


i fa ee 
lim inf = 


no log|a,, 


Quart. J. Math. Oxford (2), 7 (1956), 96-99. 
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whereas A = p = 2. It is easily seen that (5) will not hold for any series 
with gaps or latent gaps. I prove the following 


THEOREM 1. If (i) f(s) = Sa, en 
1 


is an entire function of lower order X (0 < A < a) such that 


(ii) log A,, — log Ans 


then A > lim inf An !OBAn 


nx logja,|-* 


THEOREM 2. Jf (i) f(s) = s a, en 
1 


is an entire function of lower order X (0 < A < 0) such that 


log n ..., logia,/a 
: ‘ D <o d ni™n+1 
co an (111) a 


n n+l 


(ii) limsup 


n= n 


is a non-decreasing function of n for n > No, then 
a 
A < lim inf 2S". 
n—-2 log ay 


We note that the hypotheses of Theorems 1 and 2 do not imply that 


f(s) is of regular growth. In fact we have the following theorem. 


THEOREM 3. There exists an entire function 


f(s) = > a,, ee 
log n 





for which (i) a, >, (ii) logA, ~ logA,,,, (iii) limsup = 0, 


no An 


i, SE ia al ee : : 
(iv) —— is a steadily increasing function of n, and (v) p > A. 
+17 “n 


n 


2. Proof of Theorem 1 
<1 for n > ny. Suppose first that 


Ss 
0< lim inf “"-°8“» = ¢ << @,; 
no Og |a,, sn 


Since > a,, is convergent, |a,, 


Then An logan >c—e’ forn > Nie’), 
log ay, es 


ie. la,,| > Ay we-©), 


Let e% = 2ANe-©). Ifo, < o < og4,, then 
M(c) > |a,, |e > la, |e%>» > ATAWlle—€DAnAalle—€) — DA, 


3695 .2.7 H 
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Hence, if logA,, ~ logA,,,,, we have 
loglog M(c) > (1—e")logA,, ., + loglog 2 
= (l—e’)(c—e’)(o,,,, —log 2)+loglog 2 
> (l—e")(e—e’)(o— log 2) + loglog 2, 
and so A > c, which holds when c = 0. If ¢c = o, the above argument 
shows that A = ~. 
Proof of Theorem 2. If 


_ logia,/a,., 


‘ —A,, : 
then (oc) = |a,, em for bn So < Wyss. 


Pnst ms a 


n+l 


Suppose first that A < oo. Then, from (3), 
log u(a+D+e) > log M(c) > e4-«™, 
so that log a,, +(o+D-+e)A, > eA-<e 
for all o > o, and for all n such that J, < 0+ D+e < y,,). 


If x= A,, logA,, 


then, for n > n,, 
X> 


= - m, 
log |a,, 


A,, logA,, 
(o+D+e)A,—e4-< 
> A,, log A,, 
iis A, (A—e’) HogiA,, (A—e’)}+(D+e)A,—A,, (A—e’) 
m (A—e’ )logA,, 
~ logA, —1—log(A—e’)+(D+e)(A—e’) 
Further, if %, = %,_, = ... = ¢,_, andif 1 < p< m,n—m > n, then 
we get from (6) that 











A,-p logA,,_p (A—e’)logA,,_, 
log|a,,_,|~ log A, _, — 1—log(A—e’)+(D+<«)(A—e’) 


Hence lim inf A 


nx logia, 





log A 


-1 


n 


If A = o the argument shows that c = oo. Hence the theorem is proved. 


Proof of Theorem 3. Let n, = 2, n,,, = nf (« = 1, 2, 3....), r = 1, 


Tm =m (n, << m < n?), 


m 
Nyeiy—™M 2 J 9 
rT... = Mn+1— 7 __\innwe a! (nz a m <= Nyesas Kk => s “> 3,00), 
(Mea) fA} 
- ens 


adam > a vate 


1 
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log|a,,/@,44 
Ansi—An 


which is a steadily increasing function of n. Also 


Then a, > 0 and 


= logr 


n+l? 


_logia,|~? _ logr,+...+logr, 


(in) = = 
A,, log A,, nlogn 








Hence 
(nt—n?)log(n) 


O(N,44) ™ 
. 4n$ log n,. 


1, 





(n? log n,—n?)log(nt)+ O(n? log n,) ai 


A{{n? log n,}] ~ : 
[in log n,5] n2 log nf log(n? log n,) 





It is easily seen that 
lim sup 6(n) = 2, lim inf @(n) = 1. 
n—-=x nD 
Hence f(s) is an entire function of order 1 and lower order 4. 
Finally, I wish to thank Dr. 8. M. Shah for suggesting the problem 
and for his useful guidance. 
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FINSLER-RIEMANN SYSTEMS 
By J. G. FREEMAN (Bradford) 


[Received 9 September 1955] 


1. Introduction 
It is possible, under certain conditions, to express the first fundamental 


form of a subspace of a Riemann space in the form 
ds* = g,,(x, u) dxtdz’, 


where the z’s are the coordinates of the subspace and the u’s are the 
components of a vector-field in the subspace, the g’s being homo- 
geneous of zero order in the u’s. 

This first fundamental form is similar to that for a Finsler space, 
and the properties of such a Riemann system are here obtained and 
compared with those of a Finsler space. 


2. The first fundamental form 
In a Riemann space S, of y = 2n—1 dimensions and coordinates y%, 
equations of form 


rh: ¢” +1 7” +2 


determine, for given values of the ¢t’s, a subspace S,, of n dimensions 
and coordinates 2*. (Affixes a, 8, y,... will run from 1 to v; a, 3, ¢,... 
from | to n; and A, B, C,... from n+1 to v throughout.) Variation of 
the t’s is accompanied by deformation of S,,. 

Equations (2.1), for constant values of the z’s, also define a family 
of subspaces of n—1 dimensions, in each of which the t’s may be taken 
as coordinates. These subspaces will be called the trajectories of the 
deformable subspace S.,,. 

In S, an (n—1)-parameter family of curves will be defined by n 
equations of form 

a as aM, M,...,0°-*; ps; OO... 8) (2.2) 
in which, for a given curve of the family, the A’s are constant and yu 


alone varies; these curves will be called the generators of S,,. If equations 
(2.2) are solved for the A’s and yp, equations of form 


A= Al(a,t), ..., AM =A" ez,t), p=pl(z,t) (2.3) 


are obtained. 
Quart. J. Math. Oxford (2), 7 (1956), 100-109. 
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Let u* be the S,-components of a vector-field everywhere tangential 
to the generators. The vector of this field at any point of S, will be 
called the element of support at that point. 

The u* will be proportional to the éx*/éu, and, if one of these does 
not vanish, say éx"/én + 0, we have 

ut/u” (Ga*/Op)/(ex"/Op) (a = 1,2,...,n—1). (2.4) 

Since the éx* @u are functions of the A’s, » and the t’s from (2.2), using 

(2.3) we may express them as functions of the z’s and t’s. Then @2*/éy 
will be of form f“(x,t) and equations (2.4) may be written 

ut/u" = f(x, t)/f"(z,t) (a = 1,2,....n—1). (2.5) 

If these n—1 equations are solved for the t’s, we shall in general 


obtain a result of the form 
t4 — {4(2, x), 2.6) 


in which the t’s are homogeneous of zero order in the w’s. 
If the first fundamental form of S, is 


g 2B dy*dy?, 


that of S, for given values of the t’s will be 
Ja, dx*dx”, 
where Jar = Jap Bt BE 2.7) 


and Be = dy*/da*. (2.8) 

The g, are, in general, functions of the y’s, and, through (2.1), of the 
x's and t’s; the B* are, in general, through (2.1), functions of the z’s 
and t’s. Hence, from (2.7), the g,, are also functions of the z’s and t’s; 
by using (2.6) the g,,, may finally be expressed in terms of the x’s and 
u’s, and the first fundamental form of S,, for given values of the ¢’s 


becomes 

[1] ds? = g,»(x, u)dx%dz’, (2.9) 
where the g,, are homogeneous of zero order in the components of the 
element of support. 

This first fundamental form is similar to that of a Finsler space,} and 
the deformable space S,,, together with its families of trajectories and 
generators, will be called a Finsler-Riemann system. 

The length L of the element of support is given by 


LT? = gap(x, uuu. (2.10) 


+ E. Cartan, ‘Les espaces de Finsler’ (Exposés de Géométrie, 79, 1934), 4 (3). 
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In Finsler geometry} we have the relations 

Jap = $27 L?/Ou%du?. (2.11) 
In general, the g,, of S,, will not be restricted in this way, and the system 


will be called a generalized F-R system; if, however, equations (2.11) are 
satisfied, it will be called a normal F-R system. 


3. The induced connexion 

If a vector tangential to S, has S,-components and S,-components 
X¢ and X* respectively, its absolute differential DX* induced in S,, 
when it is displaced in S,, for given constant values of the ¢t’s, is usually 
defined as the projection on S, of the vector DX*, its absolute differential 


in S,. Thus DX« = Be DX, (3.1) 
where DX* = dX*+ XPwg%, (3.2) 
a a 
wg* = ay, } (3.3) 
and Be = "Gap Be. (3.4) 
I shall adopt the same method of defining DX in cases where the t’s 
vary, in which, as the vector X° is displaced in S,,, the latter undergoes 
deformation. 
Since X¢ = Bek*, 
it follows from (3.1) and (3.2) that 
DX* = Ba BydX°+dBz X°+ BP Xw,"). 
Then, since Be By = 86, 
DX4 = dX*4+ X°w,*, 
where w,* = BS(d By + Bh wg?). 
Similarly, by definition, 
DX, = Bi DX, = Bi(dX,—Xguw,') 
= Bo Be dX,+dBe X,—Xgw,f) 
= dX,—dBy BY X,— BY ByX,w,8, 
so that, from (3.6), DX, = dX,—X),»,’. 
Now, when both the x’s and t’s vary, 
d By = Boy dx? + Bo, dt4, 
where Be, = My*/dx*be’, Be 4 = Oy*/dxdtA, 


+ Ibid. 11, VI. 
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and, from (3.3), w,P = (BY dx*+ BY, aus), B | 
S 7 
Using these expressions for dB% and w,? in (3.6) we get 


cyt = dat Be Bi,+ Bh Bt Ne }) +ava.Be By st BBB, 6 )). 
By By 


The Christoffel symbols ft 8 5 of S, and S, are connected by 


B yf \be 
Bl Bs. + BE BY,” \| —{,"). (3.9) 
y be 
The symbol f y will be defined by replacing c by A in (3.9): thus 
b 2 
5 + BP By, A \| = (,” a} (3.10) 
AB I 
We now obtain = = dat ‘ \ 4 atl . \, (3.11) 
\b ¢} \b AJ 
which reduces to “! when the ¢’s are constant. When the ¢t’s vary, 
bh C 
we have from (2.6) dt4 = pA dae+q4 du’, ‘ (3.12) 
where Pe = at4/dze, gd = atA4/due. (3.13) 
Then, from (3.11) and (3.12), 
a c G la A ns \ + cgA = 14 
w,? = da if ol | pd , | | dug? . a} (3.14) 
and we have finally, from (3.5) and (3.14), 
[2] DX4 = dX*+-X°dxT 4,4 X°duC,*,, 
oe a mf *\if © los 
anno Pye hg cj" \b “AjPe> (3.15) 
and C,*. = L * |e 


Thus DX¢ is similar in form to the absolute differential in a Finsler 

space.t The expressions ky \ vi are seen from (3.9) and (3.10) to 
\b cf’ \b A 

be functions of the x’s and ¢’s and are therefore homogeneous of zero 

order in the u’s, whilst the p4 and q4 are likewise homogeneous in the 

u’s, being of orders 0 and — 1 respectively. Hence the connexion coefficients 


+ Ibid. 5 (4). 
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T';%., Cp%, are homogeneous in the u’s, being of orders 0 and —1 respectively 
as in Finsler geometry.t 
From (3.13), since the ¢’s are of zero order in the w’s, 
wigs = 0, (3.16) 
Then, from (3.15) and (3.16), 
[3] ucC,*, = 0, (3.17) 
as for a Finsler space.t 
An alternative form for DX“ in terms of the unit vector (components /*) 
in the direction of the element of support can be obtained as in Finsler 
geometry. Since u* = La 
writing &* for Di* we have 
@* = dl¢+l'dxT,2,+1°(ltdL+ Ldl’)C,*, 
from (3.15). Using (3.17) and indicating contraction with | by suffix 0, 
we obtain &* = dl¢+dxT,c,+ LalC,*, 
= dl°(5¢+- LC 5.) +-dx°T 54, 
= dl°6¢+dxT,*., 
where 6¢ = 8¢+ LC,".. (3.18) 
Let ¢2 be the minor of 6¢ in the determinant @ = |0¢| divided by @, so 


we Og2 = 82, O2g5 = BF. (3.19) 
Then Hpt = dl°s¢+dxT,*, 4, 
i.e. dl4 = (&*—dx°T,*,)¢4, 
whence dut = Ldl¢+l4dL = L(e*—dxTy*,)¢2+14 dL. 
Substituting this expression for du“ in (3.15) we get 
[4] DX4 = dX4*+X°dxIT*,*, re 
where r*,*, = Ty*.—A,% Te’ , (3.20) 
and A,?, = LC,%, 9h | 
The first two sets of these equations hold in Finsler geometry,§ in which 
A,*, = LC,*,, corresponding to the third set.|| 

In general the [’*,*, will not be symmetrical in the lower affixes, and 


in this respect they differ from the corresponding coefficients of a 
Finsler space.tt 


+ Ibid., from 11, V and 15, (19). ¢ Ibid., from 11, V. 
§ Ibid..17 and 14 (12). | Ibid. 13, § 9. tt Ibid. 10, E. 
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From (3.17) and (3.18), uct = u*, 
From (3.19) and (3.21), ude go = ue, 
whence utg? = u?, 
From (3.17), (3.20), (3.22), 
- ; u°A,*, = ue LO,*, dh = uw LC,*; = 0. 
Therefore 
[5] A, = 0 (3.23) 
as in Finsler geometry. 

When DX¢ = 0 and S,, undergoes deformation, the vector with com- 
ponents X¢ will be said to be transported by induced parallelism in S,, 
with deformation of S,,. 


4. The generalized D-operator 
For an expression T 3b: containing affixes of S, and/or S,, we define 


DT 562: = AT 3:0: + TR ts: ey — T5.5:: wp? + --- 

+ Th 6 8 — Teper oye... (4.1) 
so that, as particular cases of this operator we have DX“, the absolute 
differential of a vector in S,, and DX¢, the induced absolute differential 
of a vector in S,, when the latter may, or may not, undergo deformation. 
With the definitidn of (4.1) the usual formula for differentiating a product 


holds, and, from (2.7), 
Day = Dag Bi By + Gag DB: By + Gag Bi DB. 
Now Dg,g = 0 and 
Jug DB3By = Jagd B3+ BY w,*— BE "|B, 
= Joy Bid BY + By w,*— BE w,') 
= Foy BEd BE + BY w,*)—w,4) 
= Inglwa!—w,!], 
from (3.6). Hence Jag DBt BP = 0 (4.2) 
and Dg» = %. (4.3) 
When g,, is expressed as a function of the x’s and t’s, I denote partial 
differentiation of g,,, with respect to 2° and t4 by (@9,,/@a*), and @9,,/ét4 
respectively. When the ¢t’s are expressed in terms of the z’s and u’s 
through equations (2.6), so that g,, becomes a function of the z’s and 
u's as in (2.9), I denote partial differentiation of g,, with respect to 2° 
+ Ibid., from 11, V. 
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and u® by (ég,,/@2*), and ég,,/@u° respectively. Then (4.3) may be 


written 
49 ap —Gara'—IJaay? = 9, 


i.e. d2*(€9 ap, C2“), + dU" (Og ap/Cue)— 
—Jap(dxT,4,+-duC,*.)—JaaldxT,4,+-duC,%,) = 0 
for all values of dz and du*. Hence 
[6] (2G ap/@X*)y, = Tove t+Uyac | (4.4) 
Jan ou! = = Crr. ae iF 
as in Finsler geometry.* 
If X is the length of vector with S,-components X¢, 
(X)? = J. XexX 
and d(X)? = D(X)? = g,, DX* X°+9,, X* DX° 
from (4.3). Hence, if DX¢ = 0, d(X)? = 0, and the length of a vector 
remains constant if the vector is transported by induced parallelism in S,, 


with deformation of S8.,,. 
Further, if X*, Y* are the S,-components of two unit vectors and @ 


is the angle between them, 
cos 6 = g,, X*Y° 
and d(cos @) = D(cos@) = g,, DX*Y°+g9,, X*DY°. 
Hence, if DX* = 0 and DY® = 0, d(cos@) = 0 and the angle between 
two vectors tig? constant if the vectors are transported by induced 


parallelism in S,, with deformation of S,,. 
This, and the previous property are posnened by vectors transported 


by parallelism in a Finsler space.{ 


5. Orthogonal trajectories 
If Jug Bz BE, = 0, (5.1) 
the trajectories are orthogonal to S,. Differentiating (5.1) with respect 
to x gives 
(29,.9/@x”) By Bz BE. +9, Bi, Bo, +9ag Bz Be. = 0. 
Writing the first term as 


([xy, B)+[By, «]) BY Bz BS, = ({ SP hate | {2° }da0) Bt Bs Be 


=| * \7,, By Be B3+| P Jose Bt Bs Bt 
p 


yf}? ‘|p 


+ Ibid. 5 (7). ¢ Ibid. 5, § 2. 
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we now have 


gop B4( Bi T BE B| : }) + 9a Be( Boot B 
\ PY 


in which the second term is 


Jan Bs{ BA. + Be, Bt By! {P= = ta a 


from (3.10). Then 


et acis Ba (Bs: Be By * ) 
Jad). Al Jap Pa act DE "\p yl . 


From (3.15), 


y = 
bde 9 ac a| 


4 \9 oa il B + Be Bz! * | 
4} tf = —Iap BA Bia BS Bi, Na 


which is symmetrical in b and d. Hence, if the trajectories are orthogonal 
Crae _ Care — $694» uf 
by (4.4); these equations hold always in Finsler geometry. 
If P, Q, R are neighbouring points of S,, with coordinates x*, 2*+-52%, 


+Az*, the S,-components of vectors PQ and PR are 527 (= X* say) 
rt Az* (= y “ say). If, as S, undergoes deformation, the points 
P, Q, R travel aldng trajectories, their S,-coordinates remain constant 


—> _> 
and the components X* and Y¢ of vectors PQ and PR also remain 
constant; the vectors are then said to be displaced uniformly. For such 
displacement we have 

Jad X*DY® = Jad X*Y dul,’ ,, 
i.e. X,DY° = X*Ydu%C,,4 
Hence, if Crue = Corer X,DY* = Y, DX°. 
Thus if the trajectories are orthogonal 


—_> —> 


X.DY =Y.DX (5.5) 


, : J om 9 
for uniform displacement of X and Y. 


A similar formula obtains in Finsler geometry{ for two vectors when 
the element of support varies whilst the components of the vectors and 
the coordinates of their common centre remain fixed. 


+ Ibid. 11, IV, V. t Ibid. 10, C. 
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6. System with deformation vectors parallel in S, for displace- 
ment along generators 
The n—1 vectors with S,-components BY we call the deformation 
vectors. For displacement tangential to S,, 
pm anhary| =(Pev) 
and for this to vanish for displacement along a generator 


(Bia+ BE BY ele a 


But then uC, = wBt| Bu, +B? omy \ ad iat 
from (6.1). Hence, if deformation vectors are parallel in S, for displace- 
ment along generators, then wC,?, = (6.2) 


These equations hold always in Finsler hehe 
Also, in this case, from (3.18), (3.19), and (6.2), 
= 3s, of =F (6.3) 
and, from (3.20), it follows that, if deformation vectors are parallel in S, 
for displacement along generators, then 
Aya = IC,*:- (6.4) 
These equations hold always in Finsler geometry. 


7. Normal F-R systems 
If the system satisfies the conditions of both § 5 and § 6, we have 
u“Cia, = 0 
from (6.2), and Cape = Chae = $69 qp/Ou* 
from (5.2). Then, from L* = g,, uu’, 


we have 


éL*/éue = 2C,,, utu’+g,, u’+g,,Uut = 29, u", 
62 L?/6uceut = 404? + 2Gcq = 22ea- 


Hence, if (i) the deformation vectors are parallel in S, for displacement 
along generators of S,,, and (ii) the trajectories are orthogonal to S,,, then 
the F-R system is a normal one. 

The following will also be proved: given an arbitrary function 2 F(z, u) 
homogeneous of second order in the u’s, there is a normal F-R system in 
which the square of the length of the element of support is equal to this function. 


Define g,,(z,u) by Jap = &F /du%du?. (7.1) 
+ Ibid., from 11, V, VI. t Ibid. 13, § 9. 
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Define generators in S,, by equations of form 
al = e+ ter), ..., amb = fm-lyt em), 2” = P+ EA), (7.2) 
where # is an arbitrary function of the A’s and ». Then 


2 
o?/es?/.../a® = acl =e GM +h gne8/./1, (7.3) 


Through (7.3), g,,, a8 preva a Pat (7.1), can be expressed in 
terms of the x’s and t's. Replacing t4 by y4 let the expression which 
is obtained be denoted by g’,,. If we now take for the first fundamental 


form OF Sy ds? = gy dy*dy? + 2gandy*dy®+9'sn dydy®, 
where the g),, and g’,, are arbitrary functions of the y’s, and define S,, by 
y=2, yt=, 

the first fundamental form of S, (¢4 constant) is 

ds? = g,,,dx%dz", 
which possesses the property g,, = }¢*L*/éu%éu". 
8. Systems in which v = 3 

It is convenient to drop the affix notation in this case and use 
X,Y,Z_ for coordinates in S,, 
x,y for coordinates in 8,, 
u,v for S,-components of the element of support, 
t for parameter producing deformation of S,. 
If x, y, t, w are eliminated from the equations 
X = X(z, y, t), Y = Y(z,y,t), Z = Z(x,y,t), 
which correspond to (2.1), and 
x = x(A,p,t), y = y(A,p, t), 

which correspond to (2.2), an equation of form 

T(X,Y,Z,A)=0 (8.1) 
results. Generators could therefore be regarded as the curves of inter- 
section of S, with the family of surfaces given by (8.1). On a generator, 
X, Y, Z of equation (8.1) are functions of zx, y, t, and along a generator 

(eT /éx) dx + (@T /éy)dy = 0, 
i.e. (éT /éx)u+(éT /ey)v = 0. (8.2) 
Elimination of A between (8.1) and (8.2) gives an equation of form 
u/v = a function of x, y, t, 

whence t = a function of x, y, u/v 


corresponding to (2.6). 





ON THE REDUCIBILITY OF POLYNOMIALS 
OVER A FINITE FIELD 


By STEFAN SCHWARZ (Bratislava) 


[Received 18 October 1955] 


Let f(x) be a polynomial of degree m over a finite field 7 = GF(q), 
where g = p*, pis a prime, s > 1. Let 

fla) = flay fala)...f,(ary (ky > 1) (1) 
be its factorization into different irreducible polynomials over 7'. Let 
the degree of f(x) be m; > 1 (t = 1, 2,..., r). 

Recently M. C. R. Butler (1) has given in this journal a non-tentative 
method of determining the number r of different irreducible factors of 
f(x). His result is in some sense very close to the results obtained by 
K. Petr (2) and by the author (3). 

Under some supplementary suppositions the results of (2) and (3) 
make it possible to determine uniquely not only the number r but also 
the numbers m,, mg,..., m,. 

Construct the following expressions, to modulus f(z), 

2D? = Cog tq F+.--- +g. n—1 2 
8 = Cyg Cy T+... +01 7" 


— = wie m—1 
r= C20 Coy Le Cg m1 X 


BODE = Cy tga Tt HOt 2 

where c;, are polynomials in the coefficients of the polynomial f(x): 
clearly Cog = 1, Co; = 0 for i = 1, 2,..., m—1. Denote the matrix [c;,| 
by C and the unit matrix of order m by E,,. 

Butler (1) proved that the rank h, of the matrix C—E,, is related to 
the number r by the relation r = m—A/,. 

Petr (2) and the author (3) proved that the characteristic polynomial 
of the matrix C can be written in the form 

C—AE,,| = (—1)™(X™— 1)(N2— 1)... (AM — 1) Aer Dom Hea Dita to Hey Dy, 

(3) 

Petr proved this formula in the case k, = k, = ... = k, = 1. The form 
(3) is given in (3). Let us remark that in (2) and (3) only the ground 
field 7’ = GF(p) is considered. The generalization to any finite field is 
ob ious. 
Quart. J. Math. Oxford (2), 7 (1956), 110-24. 
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To every irreducible factor f,(x) of degree m; > 1 there corresponds 
in (3) a factor X%—1. Conversely to every factor A*—1 in the decom- 
position of |C—AE,,| in the form +? JT] (A*—1) there corresponds an 
irreducible factor of degree s; of (1) with the exception of two cases: 
(i) if at least one of the numbers m, is divisible by p, (ii) if at least p of the 
numbers m,,..., m, are equal. For, if in the first case, say, m, = pm}, 
then ym —1 = (Am— 1), 
If in the second case, say, m, = ... = m,, then 

(™—1)...(A"»—1) = (AP™—1). 

These exceptions cannot occur if p > m. Thus in this case the decom- 
position of the form 


C—AE,,| = T] (*—1) 


determines the numbers m,,...,m, uniquely. In general the decomposition 
of the above form determines the numbers m, uniquely if there is no 
identity in 7' of the form 


(A™ — 1 )...(A%— 1) = (Ari— 1)...(A""—1), 


where the set of integers {m,,..., m,} is different from the set {m)},..., m)}. 

In this paper I give first a new proof of Petr’s theorem. The proof is 
so arranged that it gives at the same time a proof of Butler’s theorem. 
I show also how analogous new results can be obtained (see Theorems 3 
and 4 below). Then an entirely new result of essential character is proved 
in Theorem 5. Denote by o; the number of different irreducible factors 
of degree i of the polynomial f(x) (0 < io; < m). Theorem 5 gives a non- 
tentative method of determining the numbers ay, o9,..., ¢,,, or what is 
almost the same the numbers m,, mg,..., m,. In contradistinction to the 
formula (3) this theorem gives the solution of the problem without any 
exception and does not use the decomposition of the characteristic 
polynomial. 

It should be remarked that for the case p > n, k, = ... = k, = 1 the 
author used in (4) Petr’s theorem to find explicit formulae for the 


numbers o; (¢ = 


1. In this section I quote some elementary facts from the theory of 
rings which are presupposed. 

Suppose that (1) holds. It is well known that the set of all remainder 
classes of the ring 7[x] modulo f(x) forms a ring T[x]/{f(x)}. Each 
remainder class can be uniquely represented by a polynomial of degree 
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not exceeding m—1 with coefficients in 7’. Denote the set of these 
polynomial representatives by R. Then 
R = T[x}/{f(zx)} 
and the calculations in R will be carried out in the known manner 
modulo f(x). R is the direct sum of r primary rings R,, 
R= Rk, @OR,@...@ R,, 
where R, = Tla\i{fd(x*} (= 1 
Denoting the unit element of R; by »,(x) € R, we have in R 
Mm(X)+.-.+7(x) = 1, 

nil? 4d 

(ndz) (i= 4). 
Further, in the sense of set-theory, R; = Rn;, ie 

R = Rn, @ Rn, @ ... @ Rn. 


ee 


n(X)n;(x) 


With a suitable numbering of the »,(x) we have 
_ {1 (modf,(x)*) 
Lo (modfj(x)4) (j +t) 


n(x) = 


for every + = 1, 2...., r. 
Every element p(x) ¢ R can be uniquely written in the form 
P(X) = py(X)qy (2) +... + Py(%) (2) 
with suitably chosen p,(2),..., p(x) € R. 

Consider now the ring 7[x}/{f(x)*}. Write x«;=m;k; Each 
remainder class can be again uniquely represented by a polynomial of 
degree not exceeding «;— 1. Denote the set of polynomial representatives 
of T[x]/{f<x)*} by K;, where now all calculations in K; will be carried 
out modulo f(x)*‘. Then 

Ry = Trish} = Ky 
The isomorphism of #; to K; can be defined as follows. Let p(x)n;(x) be 
a polynomial € R;. Let p(x) be the polynomial € K; such that 
p(x)n(x) = P(x) (mod f,(x)*). 

With respect to (4) this means simply 

p(x) = p(x) (mod f¥). 
Then P(x)n (x) € R; <> p(x) € K; (5) 
is an explicit description of the isomorphism R; ~ K;. This isomorphism 
has the following property: if a is any element of the underlying field 7 
(a€ 7), then an,e Rw aceT cK, (6) 
This will be used in what follows. 
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2. We may regard the ring F as a finite modul over 7 of dimension m 


with the basis : a 
1, 2, 2*,..., 2". (7) 


Let a? 2 (7a) 
be any other basis of R over 7. Introduce the column matrices 
X = col(1, z,..., 2*-), Y = col(y;, Yo,---5 Ym)- 
To the change of basis there corresponds a regular linear transformation 
with coefficients in 7’ and matrix P such that 
Y = P&. (8) 

Consider now the mapping S = {z > 2%}, which carries every element 
of & into another element of R. The mapping S is a linear transformation 
of # into itself leaving the elements of 7' unchanged: for (i) if z, € R, 
z, € R, and z, — 2, z, > 24, then (z,+-2,) > (z, +22)" = 24+-2, (ii) ifa e T, 

R, then az —> (az)? = a%z4 = az’, 
I introduce the one-column matrices 
XS = col(1, 24,..., 2-9), Y5 = col(yf, 9f,..., y4,) 
and wish to find the matrices that correspond to the linear transforma- 
tion S in the bases (7) and (7a). In the basis (7) we get simply 
. XS = CX, (9) 
where C is the matrix introduced in (2). In the basis (7a) we get 
analogously ys : C,Y (10) 
with a matrix C,. It follows from (8) and (10) that (PX)S = C, PX. 
Since the elements of the matrix P are in 7’, it is further true that 
(PX) = PX*%. Hence 
Px? = C, Pa, a* = P“C, PX 
and comparing this with (9) we gett C = P-'C, P. Therefore 
C—AE,,| = |C,—2AE,,}. (11) 
This shows that the characteristic polynomial belonging to the mapping S 
is independent of the choice of the basis of the modul R over T. 

We find therefore a new convenient basis of the modul R over T. 
The ring R; = Rn; ~ T[x]/{f(x)*} itself can be regarded as a finite 
modul over 7’, of dimension «;, where we notice explicitly that, for 
a +beT, an; ~ bn; and R; contains a sub-field 7'y; isomorphic to 7’. 
Let a basis of R; over T'n; be 


(1) (2) (i) 
€; Nis & Niseees €F . Hie 
I 
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Then the m = saab elements ¢ R 


? Ma» €2? Nas «+» €F” Nay - po sony Ny (12) 


clearly form a basis a R over T. 

In the basis (12) the matrix C, has a special form which we find as 
follows. Since »; is idempotent, the mapping S assigns to the element 
e) n, (1 <1 < x,) the element 

[ef? ni]? = [eP en? = [eP Ps = ve) € Ri, 
where p(x) = [e}* is an element € R. Since the transform is in R;, we 


can write 
[ef (2) ni |? = 0 é{)) nitBi é (2), Lit ++ + Bid fe? Ni» (13) 
where the § are in 7’. Therefore the matrix C, is of the ‘diagonal’ form 
C, = diag{/C®, C®,..., C], (14) 
where C® = [Bf] (1, k = 1, 2,..., «;,). 
With respect to (11) and (14) we get finally 
C—2E,,| = |C%—2E,,, | |\C®—AE,, |...\C—AE,, |, (15) 


where E; is the unit matrix of order i. 


3. It follows from the above considerations that the matrix C® corre- 
sponds to the linear on Sigee* induced by S = {z > z*} on the sub- 
space R; with the basis e$” n;,..., €}*? n;. 

Consider now the isomorphism of K; and R; defined in (5). Then the 
elements é{”,..., e € K;, where e n; = é (mod f#), form a basis of K; 
when we regard K, as a finite modul of dimension x; over 7. Since 
(6) holds, the matrix corresponding to the mapping S defined now on 
K; is identical with C®. This follows explicitly from (13). For (13) (with 
respect to (5)) implies [¢ n,]¢ <> [e?]2 and 

Bi nc = (BP ae 14) <> BP EY, 
and therefore [eM }e = Bi ef + PP e+... + Bid exo, 
In any other basis of K; over T the matrix D© belonging to this mapping 
is similar to the matrix C®, 

We choose a new basis of the modul K; over 7. It is clear that the 


x; elements 
}. 2, fs gmi-l 


fiz), xf{zx), sing x™-1f (x) (16) 
are, “rfp, 2, emia 


are linearly independent with respect to 7’ and therefore form a basis 
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of A,. The linear independence follows from the fact that the degrees 
increase. I shall establish the form of the matrix D® in this basis. 
Take any element 

afizyp (O<y< m—1;0<8 <k,—1). 
The transform 2*7f,(x)® ¢ K, can be written as a linear combination of 
the elements (16) with coefficients € 7’. Accordingly 


. mi—1 —<_ mi—1 mi—1 
xvIf (xa = F aPet FY aMPafi(x)t...t+ > afP_, af (x) 
l=0 i=0 l=0 
F , (17) 
with af%), af? 7) eT (l= 0, 1.,...,m,—1). 
Suppose now that § > 1 and take (17) to modulus f,(z)*+!. Since 
dq > 5+1, we get 
ie mi—1 5 3%! 3 5 
Y aR e+fiz) > aff? a'+...tflzyP J ofa? (mod f,(x)°+?). 
l=0 l=0 i=0 
Hence all coefficients nf), (7°) ax) (0 <1 < m,—1) must be zero. 
The relation (17) is necessarily of the form 
mi—1 


" mi—1 
xv4f (ar) >, a8) aif (x)P+4+...+ on al) alf(x)e 2. (18) 


The relation (18) implies that the matrix D® is of the form 


- (i) (i) 5 
Di | DS‘ 











“sO 


in which all elements in the (m;+1)th, (m;+2)th,..., (m;&,)th rows in 
the main diagonal and to its left are zeros. 
Since C™ is similar to D®, we have 


C®—E,,| = |D©—2E,, 
and, with respect to (19), 
D®—AE,,| = |D—AE,,,|(—A)™*». 
The relation (17) with 6 = 0 is of the form 


mi—1 


; mi—1 mi—1 
SBP J aM aplayt+.+ 2 abee-aefdee, (22) 


and therefore D{ = [a9] (1, y = 0, 1,..., mj—1). 
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The relation (22) shows that the matrix Df? is identical with the 
matrix obtained by calculating 
0 - . -mi—1 
a Boot+Bor *+---+Bom —1 2" 


xt = By tBu t+ --+Bi mre (modf,{zx)), (23) 


{mi—1 Q ae?) eS my -mi-l 
wt Pmg 1,0 Pm; 11% ad -* Lm-1% ; 


where I have written of,” = £B,,. 


4. The problem of finding the characteristic polynomial of |C—AE,, 
reduces to the case of finding the characteristic polynomial of the linear 
substitution (23), where f,(z) is an irreducible polynomial of degree m; 
over 7’. 

For convenience let us write (in this § 4) d(x) instead of f,(7) and n 
instead of m;. Denote by j a root of d(x) = 0 and consider the extension 


field K = T[j] ~ K[x]/{d(x)}. 

We again regard K as a finite modul of dimension x over 7’ with the 
basis fo Pyenes fh. (24) 
The mapping 


is now an automorphism which transforms the elements (24) into the 


elements 


These elements are linearly independent with respect to 7’ and there- 
fore form a basis of A. For, if there were n numbers yp, 7;,..., ¥,-1 € T; 
not all zero, such that yy+y,j?+...+-y,_,j"™ = 0, we should have 
a wey 2 ee in—1]a — ee OS aaae n-1 ' : 
[YotVJ t+ +Y¥naJ” lt = 0, Le. yotyjt+--+ynaj"* = 0, and j 
would satisfy an equation of degree less than n, contrary to the assump- 
tion of the irreducibility of 4(). 
Writing 
J = col(l, j,..., 9*~*), F = col(1, #.,..., #*-™), 
we get J* — BJ, (26) 
where the matrix B, identical with the matrix D™, is the matrix of the 
substitution (23) and is, of course, regular. 
. ° P oh oe ° k 
Denote by J* the one-column matrix withelements 1, j“, j,...,j®-2 
‘Since, for every a € T, a4 = a, we get the system of equations 


J* = BJS, J* = BI", ..., J* = BI". 


Since j satisfies an irreducible equation of degree n, we have j” = j 
«_k —* . 
and, for every k <n, jf #j. Hence JS* = J and, for every k < n, 














ON THE REDUCIBILITY OF POLYNOMIALS 117 


J** + J. The last of the equatioris (27) can be written in the form 
J = BJ*"”’. This implies successively 
I= B= BD — ... = BF = BJ. 
Hence B” E... 
The matrix B satisfies the equation (A"—1)E, = O,, 
the zero matrix of order n. I shall show that A”—1 is also the minimum 


polynomial of B and hence, with respect to the order of B, at the same 
time the characteristic polynorial of B. To this end it is sufficient to 


where O,, is 


show that the matrices 
z., B, »,.... Be (28) 
are linearly independent. From (26) and (27) there follows 
J=aE,J, F=BJ, P= BY, ..., O° = BJ. (29) 
If the matrices (28) were linearly dependent, there would exist n numbers 
,, not all zero, such that 


E,.+ 7, B+...4 Yn-1 B"-!=O 


‘Uv n 


Ver Yir-9 Vp 
= 

Multiply this relation with J and use (29). We get 
J+y,J§+y,.J3%+...+y,_,5"" = Of, (30) 
where O}, is the one-column matrix with all its elements zero. Put 
j = X%, j§ = X%, jf = Xp,..., J” * = X,_1: all these elements are different 
from each other. The equation (30) written explicitly is of the form 


70 


We | 


Yo 1+ Y1 1+-...+ Yn-i 1 0 
YotoT¥1%1 7 +++ TYn-1Tn-1 = 0 
ap_ wh—-l1sa, jana. - n—1 ' 
YoXo M17 + TY¥n-1%n-1 0 


This relation can be satisfied if and only if the system of n homogeneous 


equations 


— ape Se # ‘ 
YoCoTV1T1 + -- TYn-1%n-1 = Y (0O<l<n 1) (31) 
has a non-trivial solution in yo, y;,..., ¥,-,- This is impossible since the 


determinant of the system (31) 


is not zero. Thus we have proved 
B—AE,,| = (—1)"(A"—1). (32) 


The results (15), (20), (21), (32) together imply the following theorem: 
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THEOREM | (Petr-Schwarz). Let (1) be the factorization of the poly- 
nomial f(x) over the finite field T into different irreducible polynomials 
f,(x),..., f(x) of degrees m,,...,m,. Let C be the matrix of the linear substi- 
tution (2). Then the relation 

ss E,, _ ( a yam at )(A™— l )...(A"%— l ABD +(ke—1)m2+-"+(k,—1)m, 


holds for the characteristic polynomial of the matrix C. 


5. We have proved that the matrix C is similar to the matrix (14), 
where each of the matrices C is similar to a matrix D® of the form 
(19). Denoting the similarity by the sign ~ we have 

C ~ C, = diag/C™, C®,..., C”] ~ diag[D®, D®,..., D™]. (33) 

The minimum and characteristic polynomials of the matrix D{® are 
identical as far as the factor (—1)”* and equal A"*—1. The invariant 
factors of D{?—2E,,,, are therefore X"*—1, 1, 1 (m;—1 times). Con- 


sider the m; x m; matrix 


0 


The invariant factors of N{?—AE,,, are clearly again \”"‘—1, 1, 1 


(m;—1 times). Hence the matrices D{? and N{ are similar. There is 


therefore a regular matrix Q{? such that 
@M HMWOrQg@)] — NW 
Q;’ DY[Qy)}* = Ny. 
Transforming the matrix D® by the regular matrix 
(i) 0 
QW . Q) = —| 
0 Ein dics DV 
we get a new matrix, which will be denoted by N®, 
, N(? Ni? 4 


QOD of Q”} S st 








In the matrix N® all elements in the (m;+1)th, (m;+-2)th (m,; k;)th 
rows in the main diagonal and to its left are zeros. The result (35) implies 


diag[D®, D®..., D®] ~ diagfN®, N®,..., N]. (36) 
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The formulae (33) and (36) imply the following statement, which is the 
only result needed in the following: 


The matrix C from the relation (2) is similar to the matrix 
N = diag| N®, N®,..., N®], (37) 
where each of the matrices N® is of the form given in (35). 


The matrices C—AE,, and N—AE,, are similar: in particular they 
have the same rank for every A. 
(A) Choose A = 1. The matrix C—AE,, is of the form [see (2)] 


0 0 0 . “ ” 0 
Cio 1 ] Cie 0 (38) 
Cm-1,0 Cm-1,1 Cm-1,2 Cm -1,m-1— l 


The matrix N—AE,, takes the form 


N—E,, = diagjN®—E,,, N°—E,,,..., N°—E, ], 





where Pi EES re 0) | 
O=—332 wo sc a OG 
C Ni 
6 #26. <« <= i : 
; 1 oo ...0—-ij 
N® E, — — | 
i —l ., 
‘Gage 
0 % 
os 
= =" 








The rank of N—E,, equals the number of linearly independent columns 
of N—E,,. A column that intersects the block NO—E,, is clearly 
linearly independent of the columns that intersect the blocks NY)—E,, 
when i #j. Hence it is necessary only to establish the number of 
linearly independent columns in NO—E,,.. The (m,;+-1)th, (m,+-2)th.,..., 
(m,;k;)th columns are linearly independent and independent of the first 
m, columns in NO—E,.. It remains to establish the number of indepen- 
dent columns in N{?—E,,,. But this matrix has clearly the rank m,;—1. 
Therefore NO—E,,., is of the rank m;k;—1. Thus the rank h, of the 
matrix N—E,, is 


hy = (m,k,—1)+(m, k,—1)+...+(m,k,—1), 


i.e. hh, = m—r. We have just proved 
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THEOREM 2 (Butler). The rank h, of the matrix (38) is equal to m—r, 
where r is the number of different irreducible factors of the polynomial (1) in 
T and m the degree of f(x). 


(B) Choose A = —1. Then the determinant of the matrix 
5 = oe ° - & 8 
2 ee eS ce an eee 
NiO. E,.. _ ak a Ser es (39) 
0 0 0 ak ; ; 
2 oe es ae 








is 1+(—1)"*—. Suppose p > 2. Then, if m; is odd, the rank of this 
matrix is m,; if m,; is even, the rank of (39) is m;—1. Therefore the 
rank h’ of the matrix N+E,, is given by 


h' = >’ mk, + >” (m,k;—1), 


where >’ runs through those summands in which m; is odd and >” 
through those summands in which m; is even. We have proved 


THEOREM 3. Let p > 2. Let r’ denote the number of different irreducible 


factors of even degree of the polynomial f(x). Then r' = m—h’, where h’ 
is the rank of the matrix 


2 0 0 ae ae 0 
C10 Cy,+1 Cie a C1.m-1 
Cm-1,.0 ©m-11 ©m-1,2 © = © Qype ti 


Remark. If p = 2, the argument just used leads again to Theorem 3 
and does not give a new result. 

(C) The result just given can easily be generalized. The matrix 
Nj? —AE,,,, has clearly the rank m;— 1 or m,; according as the determinant 
N{?—AE,,,,) = (—1)™(A"*—1) is or is not zero. Let now g be any 


element of 7. Suppose that g belongs to the exponent /, i.e. g = 1, in T 
but g = 1 for every 0 <I’ <1. Then it is well known that 


g-1=0 (0<t<q—1) 


if and only if J\t. Therefore the matrix N°—gE,, has the rank «x; 
except where /|m,. In this latter case the rank of N°—gE,, is x;—1. If 
there are r” different irreducible factors of degrees divisible by /, then 
the rank of N—gE,,, and hence of C—gE,,, is m—r”. This proves the 
following theorem. 
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THEOREM 4. Let g be an element of T belonging to the exponent I 


(0 <l <q—1). Let the set {m,, mg,...,m,} have the meaning defined by (1). 
Let r° be the number of integers of the set {m,, mg,...,m,} which are divisible 
by l. Then r” = m—h", where h” is the rank of the matrix 
l—g Vv ee 0 
Co “-9 + =: - Cim-1 
Cm-1,0 ©m-1,1 . . . Cm—1,m-1— 9 


6. The result of Theorem 2 can be reformulated as follows. Let a; 
denote the number of different irreducible factors of the polynomial f(z) 


of degree j. Let the rank of C—E,, be h,. Then 
h, = m—(o,+0,+...+0,). 

In this formulation the result can be generalized in such a manner that 
it leads to the very remarkable result given in Theorem 5 below. 

Suppose that f(x) have an irreducible factor of degree j, i.e. ¢; > 0. 
In this case the matrix (37) contains a block N®, where N{° is a matrix 
of the type (34) having exactly m; = j rows and columns. 

Consider the /th power C! of the matrix C (J >1). Then C’~ N’ 
In this last matrix one of the blocks is the matrix [N} of order 





m, k jk; The matrix [N®} is of the form 
. [N{?F M, |] 
ae Se 4 
[NOF = 0 p> (40) 
0 | 
| @o-« 
where | N{?} is a square matrix of order j of the type 
(+1) 
eh eee se eee 
owe 2 | 2s ae SD 4 ee 
0 0 Se Ow « «ww S 
NO} b <abe oe 41 
IND] ei a a ae oe ee (41) 
SS. Sno eee a, Se, 
ioe . it ye eee. se 








Since [N{”}) = E; we have 


[NPP = [NP 


if p = o (mod )). 








122 STEFAN SCHWARZ 


Let us ask what is the rank of N'—AE,, for A = 1. Since 
N! = diag{[N®}, [N@F}...., [NF 
and [N®} is of the type (40) it is sufficient, as above, to consider the 
The relation (41) implies 


square matrix [Nj }}—E,. 


(l+1) . 

i" en l Bs ok ) Lee 

Swe @ es 2 0 Sop ait ag oe ee 

— 0 0 at 0 0 0 
N®}?_E. = a ies 9 ee . (42 
rT—s; , Oe. ss or el. Le (*°) 

0 l 0 QO —] 0 

. § on 0 0 .. kl 








I shall show that the rank of the matrix (42) is j—(/,j), where (/,)) 
denotes the greatest common divisor of / and j. 
Let us introduce the indeterminates z,, z,, Z3,... with the convention 


that z; = z, if i = k(modj). Then the rank of the matrix (42) is the 
same as the number of linearly independent forms in the system 
2g t-294:39 — Ze tpg, Ze 247490-°%s 2524; (43) 


To find this number we proceed as follows. Take first the forms 


~ 


“3J+1°"*** 


“1 Asa» 41 Tats F204 
In this sequence we return to the first form after « steps, where « is the 
least positive solution of the congruence 


«kl+-1 = 1 (mod/)). 
Hence « = j(/,j)-'. Then the « forms 


“1 T “1+1> 244 Zo741 geees — 2 ye 94 T Mad +1 (44) 
have the property that the first x—1 are linearly independent whereas 
the last is their negative sum. If « <j, we take the forms 


21 1+2> —“t+2 


Z Fetinseess Ze —1848 Td 48° (45) 
Here again we have « forms of which the first «—1 are linearly inde- 
pendent, whereas the last is the negative sum of the others. Hence it 
is clear that the «—1 forms from (44) are linearly independent from the 
forms in (45) since in both sets only different indeterminates z; occur. 
Repeating this process we find that the forms (43) can be divided into 
jj(l,j)* = (Lj) subsets in each of which just one form is linearly 
dependent of the others and the forms in different sets are linearly 


independent. To prove this last statement it is sufficient to show that 


vittr=pl+o (mod)) (45 a) 
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is impossible with r + a (0 <7 < (l,j),0 <o <(I,j). In fact (45a) would 
imply vlit+r =ypl+o (mod (I,j)), 

i.e. r = a (mod(/,j)) and therefore s = o contrary to the supposition. 
Hence (43) contains j—(/,j) linearly independent forms and therefore 
the rank of the matrix (42) is j—(I,j). 

If in the factorization of f(x) there is a factor of degree j, i.e. ¢; > 0, 
of multiplicity £;, then there is in N'—E,, a block of order jk; having 


the form [NO P_—E, M, 





[NOF-—E,, = 0 ae 


Accordingly the rank of this matrix is k;j—(l,j). 

Denote by h, the rank of the matrix C’—E,,, which is the same as the 
rank of N'—E,,. If there are o; different irreducible factors of degree j 
with multiplicities k}”, k}*,..., ki”, then the contribution of these factors 
to the rank is ‘ ‘a 

EPID} =I TEP—Cidoy 


the sum being empty if ¢; = 0. Therefore 


7 4 Co | m oc | m 
h, > ijI> (1, j)o;; = y i J > &Y) —¥ (,j)e;. 
1 1 j=l1' fs j=1 
m Cj 
Since > >} ki) = mM, 
j= s=1 
m_ F 
we get > (L.j)o; = m—h,. 
j=1 


Writing this equation for / = 1, 2,..., m we get the following system of 
linear equations 


(1, ljo, +(1, 2)o,+(1, 3)o,+...+(1, m)o,, = m—h,, 
(2, 1)o, +(2, 2)o,+(2, 3)o,+...+(2, m)o,, = m—hyg, 
(3, 1)o, +(3, 2)o,+(3, 3)o,+...+(3, m)o,, = m—hs, 
(m, 1)o,+(m, 2)o,+(m, 3)o,+...+(m,m)o,, = m—h,,. (46) 


This system has always a solution since it is well known that the deter- 


minant ti es ere 
Sh @h.... &a 


(m,1). (m,2) . . «~ (o,m) 
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is different from zero, and in fact equal to 4(1)4(2)...d(m), where (2) is 
the Euler function. We have thus proved 


THEOREM 5. Let f(x) be a polynomial over the finite field T. Let a; 
(1 <j <m) denote the number of different irreducible factors of degree j 
of the polynomial f(x). Let C denote the matrix of the linear substitution 
(2) and h, (1 <1 < m) the rank of the matrix C'—E,,. Then the numbers 
1, Fg,.+-, O, are uniquely determined by the system of linear equations (46). 


m* 


Theorem 5 enables us to find from a unified point of view a number 
of more or less known results concerning the reducibility of some special 
congruences. The details can be left to the reader. 
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A NOTE ON UNIFORMLY DISTRIBUTED 
SEQUENCES 


By P. B. KENNEDY (Cork) 
[Received 12 November 1955] 


1. Ler {x,,} be a sequence of numbers lying in the interval (0,1), and, 
for0 <a <8 < 1, let N(a,8;n) be the number of z, (1 < r < n) lying 
,B). If at 
whee N(a,B;n)~ n(B—a) (n>) 

for every (a, 8), then the sequence {z,,} is said to be uniformly distributed 
in (0,1). It is known [Hardy (2), 115] that, if {z,} is uniformly distri- 
buted in (0,1), then 


1 
f (2p) - > | f(x) dz (C,1) (n+o) 
0 


for every f(x) which is bounded and R-integrable in (0,1). In particular, 
if we take f(x Z, . 
air % 2, —>4 (C,1). 
On the other hand {z,} obviously diverges. Hence by a well-known 
Tauberian theorem due to Hardy [(1); (2), 121], 
lim sup n|Az,,| = +90. (1.1) 
nox 

It was shown by Littlewood [(3), Theorem (C)] that (1.1) cannot be 
improved for general divergent sequences which have (C, 1)-limits. The 
sequences constructed by Littlewood are not, however, uniformly distri- 
buted, and, as far as I can see, no slight modification will make them so. 
It may therefore be of interest to show, by a class of simple examples, 
that (1.1) cannot be improved even for uniformly distributed sequences. 
This is the content of the following theorem. 


THEOREM. Let d(n)—> 00 steadily as noo. Then there exists a 
sequence {x,}, uniformly distributed in (0,1), such that 


Uns 


nAx,, = Odf{(n)} (n>). 


+ I use the notation 
Au, Un—Un_1> Af (k) mene Sf (k)—f (k—1). 


Quart. J. Math. Oxford (2), 7 (1956), 125-7. 
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2. The proof of the theorem depends on the following lemma. 


Lemma. Let ¢(n) be as in the theorem. Then there exists a sequence of 
integers {p,}, tending steadily to infinity as k -> co, such that 


Po = 9, Ap, > ©, =. >I, ~ = Ofb(p,)}- 

Let Y(k) = max[2, min{kt, d(exp k*)}], (2.1) 
so that ¥(k) > c steadily as k > oo. The first step in the proof of the 
lemma is to construct u(t), tending to infinity steadily, so that 

wk) < o(k) 

and Apu(k) = O(k-*). 
Choose an increasing sequence of integers {m,} such that 

mM, = 0, m, = O(Am,) (ro), (2.4) 
and %(m,.) >r (r>1). (2.5) 
This can be done by making m, tend to infinity sufficiently rapidly. 
Define p(k) to be r when k = m,, and by linear interpolation for all 
other integer values of k > 1. Then, when m,_, < k < m,, we have, 

ro p(k) <r < Y(m, 4) < wh), 

which is (2.2). Further, when m,_, < k < m,, by the definition of y:(k), 





_ Kulm, )—p(m,_ i} — 
kAp(t) Am, <a 


which gives (2.3) if we make use of (2.4). 
Next we define a sequence {e,} which, except that it is not a sequence 
of integers, is essentially the sequence {p,} of the lemma. For k > 1 put 


k 
&. = exp(75}- (2.6) 


By (2.1) and (2.2), e, > exp(k!) (k > 4). 2.7) 
Further, by a simple calculation involving the use of (2.3), we obtain 
from (2.6) ey ~w pli) Ate. (2.8) 
It now follows from (2.1), (2.2), (2.7), (2.8) that 

A ose 

€k-1 


ai = OfMh)} = Ofpfexp(k4)}] = Of6(e,)}. 
k 


The conclusions of the lemma follow readily from this if we put p, = 0 
and define p, (k > 1) to be the least integer not less than e,. 


Ae, > 0, 
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3. We can now define the sequence whose existence is asserted by 
the theorem. Throughout this section the integers k and n satisfy 
2p. <M S 2Pye 13 (3.1) 
it follows from the lemma that this implies 
ui~ 2p}. ~ 2Pp41- (3.2) 
For all k > 0, put 


n—2p, 
~—“Pe (2p, << 2 < P_t+Pesr), 


Ap,. 
x, weve (3.3) 
2P eu e 
—. (Pe+Prsr < % S 2p). 
Obviously 0 < zx, < 1 for all n. 


To show that {z,} is uniformly distributed, denote by N;,(«,8) the 


number of z,, satisfying (3.1) and lying in (a,8). From (3.3) it follows 
without difficulty that 
N,.( x, B)—2(B—a)Ap,., < 2, 
whence, since Ap,., > ®, 
Na, B) ~ 2(B—a)Apy +1. 
Therefore, by a familiar theorem on limits, 


“ k-1 
N (a, B; 2p,) = 2 NilaB) ~ 2(B—a)p,, 


and so, for any n, 

N(a,B8;n) = N(a,B; 2p,)+O(Ap, 41) ~ m(B—a); 
here we have used (3.2) and the fact that Ap,,., = 0(p,), which is a simple 
consequence of the third statement of the lemma. Hence the sequence 
{x,,} is uniformly distributed. 


n) 


Further by (3.1), (3.2), (3.3), and the lemma, 





s n 2Pe+t 
n|Az,| = < = Ofd(p,.,)} = Ofd(n)}. 
in.” Bie {$(Pe+1)} {d(n)} 


This proves the theorem. 
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DECOMPOSITION SPACES AND SEPARATION 
PROPERTIES 


By M. SHIMRAT (Jerusalem) 
[Received 10 January 1956] 


Ir a topological space S is decomposed into disjoint sets S, (a ¢ A), we 
have a ‘decomposition space’ S* defined as follows. The points of S* 
are the individual sets S,, and a collection {Sg} (8 ¢ B) of such points is 
an open set in S* if and only if the union (J Sz (f € B) is open in the 
original space S. The space S* may be regarded as obtained by identi- 
fying the points of each S, into a single point. 

The decomposition space S* will in general have no separation proper- 
ties even if S is assumed to possess such properties. If, for instance, 
the cartesian plane £*? with coordinates x, y is split into the two sets 

Sy: (x,y); 2 < Of U (0,0), 
S,: £?—S,, 
then clearly the decomposition space does not have property 7). 

It is the object of this note to prove the theorem: 

THEOREM. Any topological space is homeomorphic to a decomposition 
space of some Hausdorff space. In other words: the most general topological 
space may be obtained from some Hausdorff space by a process of identi- 


fication. 


In order to prove the theorem let X be any topological space. We deter- 
mine a Hausdorff space Y which is the disjoint union Y = (J Y, (x € X), 
where each Y, is dense in Y. (To show that such Y and Y, exist, consider, 
for example, the topological product of real lines E,, one for each point x 
of X, and take as Y the set of all points of [T] £, (x €¢ X) having exactly 
one rational coordinate and as Y,, for any definite x of X, the set of 
those points of Y whose rational coordinate has index .) 

Consider now, in the topological product X x Y the set 

Z =Ute <Y,) = {exy; ye Y,}. 
Z is a Hausdorff space (in the relative topology), i.e. any two points 
of Z have disjoint neighbourhoods. This follows from the fact that the 
Y-coordinates of any two points of Z cannot coincide and that Y was 
assumed to be Hausdorff. 


Quart. J. Math. Oxford (2), 7 (1956), 128-9. 
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I shall show that the decomposition of Z into the sets Z, = xx Y, 
gives rise to a decomposition space Z* homeomorphic with X. Let us 
map every point Z, of Z* into the point z of X. This is clearly a (1, 1) 
‘onto’ mapping. If G = {Z,: 2 € G’} is any open set in Z*, then by 
definition (J Z, (x € G’) is open in Z. Let x, be any point of G’ and zy x y 
one of the points of Z,,. Since (J Z, (x € G’) isopenin Z, some neighbour- 
hood ZO (Vy(2_) x Vy-(y_)) of 2 X Yo is contained in it, where Vy(z,), Ve(Yo) 
are neighbourhoods of 2», y, in X, Y respectively. If now z’ is any 
point of Vy,(z,), then the set Z, = 2’ Y, must contain a point of 
Z 1 (Vx(xq) x Ve (yo)) since Y, is dense in Y. Hence, by the above, 
Z,, meets [J Z,(x € G’), and so the point Z, of Z* lies in G. Its map 2’ 
therefore lies in G’. We have thus shown that G’, which is the image 
of G under our mapping, is open. 

On the other hand it is obvious that, if G’ is any open set in X, then 
its counter-image {Z,; x € G’} is open in Z*. We have thus shown that 
the mapping is a homeomorphism, which completes the proof. 


3695 .2.7 EK 











A PROBABILITY THEORY OF A DAM WITH 
A CONTINUOUS RELEASE 
By P. A. P. MORAN (Canberra) 


[Received 27 January 1956] 


Suppose that a dam has capacity A and that at the beginning of each 
year ¢ an amount X, flows into it. We suppose that X, has a known 
probability distribution and that the random variables X, for different 
years are independent of each other. If Z, is the amount of water in the 
dam just before the input, the dam will contain, after the input, an 
amount X,+Z,, or AK, whichever is the lesser. An amount Y, is then 
released according to a definite rule, leaving an amount Z,,,. The main 
problem in the theory of dams is to calculate the probability distribution 
of Z,,, when the system has settled down to a stationary behaviour, so 
that this probability distribution is independent of ¢. 

To make any progress with this problem it is easiest to begin by 
assuming that X,, Z,, A, Y, take only values which are integral multiples 
of some unit. X, and Z, then have discrete distributions. A variety of 
release rules may be chosen, but the one I shall consider is that the 
amount released, after the input, is X,+ Z, or M, whichever is less, where 
M is a fixed multiple of the unit. 

If {p;} (i = 0, 1,...) is the probability distribution of X, and {P;} that of 
Z,, we then have the system of equations (5) 


Py = P(pot---+Pu)+Py(Pot---+Py-1)+---+Pay Po 

Po = hypys T41 Pu ot Py Pit Py si Po . (1) 
Prev == P(pe+...)+... the _y(Pyt---) 

In this case 0<Z,< K—M, 


and the distribution of Z, can be found by solving equations (1). This 
system of equations may be regarded as an approximation to a con- 
tinuous situation for which a corresponding integral equation can be 
found, but the latter is not easy to solve. Given any continuous proba- 
bility distribution for X, we can approximate to it by a discrete distribu- 
tion and then solve equations (1). The solution will be an approximation 
to the distribution which is the solution of the continuous problem. This 
is a practical method for a numerical approximation, but a theoretical 
solution for the continuous case has so far only been carried out when 


Quart. J. Math. Oxford (2), 7 (1956), 130-7. 
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the input has a negative exponential distribution (6) and the solution 
is then very complicated, having different analytical forms in different 
parts of the range (0, K—M). 
It would, however, be very desirable to have a solution for the case 
when X, has a type-III distribution of the form 
= : () ) A XimqX., (2) 
mI (k)\m 
The purpose of the present paper is to find such a solution in a modified 
problem where the time variable is taken as continuous and the release 
is at a constant rate. In order to do this we consider the above system 
of equations when they are modified, first by making K infinite and then 
by making VM i. 
Making K infinite we get the infinite system of equations 


Py Py Po T et py) +... +Pry Do 
P, Py Puss +...+Py py t+ Pry. Po \ (3) 


These, as will be seen, correspond to the case of an infinite dam when 
we suppose the probability density to be concentrated at the bottom. 
These are the equations of a Markov chain with an infinite number of 
states, and for them to have a meaningful solution certain conditions 


« 


are necessary. 
Now write m = > np,, which is the mean of the input. Then I shall 
0 


show that, if m < M, there is a meaningful solution, whilst, ifm > V, 
no such solution exists. The matrix of coefficients of the P; on the right- 
hand side of the above equation is, in fact, the transpose of the transition 
matrix, (p,;), of the Markov chain. We have then 


x 


> JPis -i—(M—m) (¢> MM). 
7 =C 
It then follows by a slight extension of Theorem 2 of (2) that the process 


is ergodic. 


If, on the other hand, m > M, consider 
| Zing = J t+X,—Y, 
Then Y,< ¥, E(Z,.,)—E(Z,) > m—M, 
so that E(Z,.,)—E(Z,) > E(X,+...+X,,,-,)—sM. 


Now by Khintchine’s weak law of large numbers 


a-1(X,+-...+X,45-1) 
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converges in probability to m as s increases, and so the right-hand side 
of the last inequality can be written as 3(X —M), where X —M converges 
in probability to m—M > 0. The system is therefore dissipative and 
there is no meaningful solution. This can also be made to follow from 
Theorem 4 of (2) (I am indebted to a referee for this remark). 
This type of argument is familiar in the theory of queues [see Kendall 
YI - . 1 
(3)]. We now suppose that VM = 1. Equations (3) become 
> > 2 
Fy = Fy(Pot+Pi)+F Po 
> > > > 4 
P, = ype TA PTE sPo | (*) 
These are very similar to the equations used by Kendall (3) in the 
theory of queues. The matrix of the coefficients in (4) is almost lower- 
triangular, and this enables us to use generating functions. We assume 


m= Mm, >) ip<il=WM 


0 


and write m, = > ®p; 
3 = ‘ 


for the sth moment of the distribution {p,;}.. We also write 


Plz) = > p,2', P(z) = > P, 2". 
0 0 
Multiplying equations (4) by z, z*,... and adding, we get 
zP(z) = p(z)P(z)—p)(1—z)P,: 
and so PG) = a» F, Le - (5) 
z— p(z) 


We can find P, by noticing that 


1 = P(1) = p, Pylim2—* 


2+1 2— p(2) 
= pofitl—p'(1)}, 
and so Po Py = 1—p'(1) = 1—m,. 


This result is remarkable in that it shows that the probability of the 


a 
dam being empty depends only on p, and m,. Writing VM, = > i*P, we 
0 
can show that ° 
— p'(1) 
M, = P’(1) = >——__|-°: 
21—p (1); 


and higher moments can be found in a similar way. 
We now apply these results to the continuous problem by letting the 
unit in terms of which the above discrete distributions are defined tend 
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to zero, and we suppose that the time interval between successive inputs 
also tends to zero, so that we are then concerned with a case in which 
the water is released at a constant rate. We take this rate to be unity 
and n an integer such that n-" is the unit amount released at the end of 
successive intervals of time of length n-'. The total input during an 
interval (¢,¢+1) then has a probability distribution which is the nth 
convolution of {p,;} and has the generating function {p(z)}". 

We now have to let {p;} vary as'n increases so that its nth convolution 
approaches the continuous distribution we are going to use for the con- 
tinuous input during unit time. 

Consider first the geometric distribution 


p, = (l—r)r' (¢ = 0, 1...) 
where 0 <r < 1. Then 
p(z) (1—r)(1—rz)-", p’(1) = r(i—r)-. 
Let X be a variate with this distribution. Its characteristic function is 
(1—r)(1—re*®)-!. 
Then the variate Y = n-1X has mean 
rn-'(1—r)-}. 
We let » tend to infinity and r to 1, so that this mean tends to m,, say, 


by putting 
YI 6 r= m,n(l+m,n)-}. 


Then the characteristic function of Y is 
~t 
, m,n 
¢(8) lim (m, ny] —|( mn) iain 
n—x \ l+m,n j 
(1—m,i6)-}, 
so that Y has the asymptotic distribution 


—1 
my, 


e-vim dy. 


We shall take this as the asymptotic distribution of the total input 
during any interval (f,t+1) and we do this by supposing that the nth 
,} is the above geometric distribution. Thus the distri- 
bution {p,,} must be the ‘nth-root convolution’ of the geometric distribu- 
tion and has the generating function 


convolution of { p 


p(z) = (l—r)*#"(1—rz)-™*, (6) 
where r= m,n(1+m,n)-}, m, <1. 


That this gives a genuine probability distribution can be seen by expan- 
sion. It is, in fact, the negative binomial distribution. Inserting the 
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distribution in (5) we obtain the generating function of a discrete 
distribution {P;}, which we may regard as an approximation to the 
solution of the continuous problem. In this continuous problem the 
water flows out at unit rate whenever there is any water in the dam 
and the input during unit time is distributed as 


mz te-Xim gx. 
Thus the input during any time T' is distributed as 
l X | T-1 : " _ 
en —| e-Xim dX, (7) 
m,1(T)\m, 

The total amount which has flowed into the dam from time zero to 
time 7, say, is therefore the variate of an additive random process 
which might be called the ‘gamma’ or ‘type-III’ process because the 
variate has a type-III distribution. This is a particular case of more 
general additive processes which are discussed by Lévy (4) under the 


name of ‘processes with movable discontinuities’. 
Now for (6) we have 
Po = (1—r)"™*, m, = rn—(1—r)-}, 
so that the probability of the dam being empty is 
I—p’'(1) = 1—rn-(1—r)" 
De 7 (l—r)"™ 





and, as n tends to infinity with r = m,n(1+m,n)-", this tends to 1—m,. 
We therefore suspect that 1—m, is the probability that the dam is 
empty in the continuous case. This is in fact true, as we shall see, but 
does not follow from the above since the limiting concentration of 
probability at the origin in the continuous case might have arisen in 
part from banking up, in the limiting process, of the contributions of 
P,, P,, and other values of P; with ¢ small. 

We now find the characteristic function of the distribution of Z,. 
Notice that we have defined the solution in the continuous case to be 
the limit of the discrete solution and not as itself the solution of a 
functional relation. The establishment of some kind of integral or 
integro-differential equation for the solution in the continuous case 
would be interesting, but I have not so far been able to do this. 

Writing y = e’? we obtain for the characteristic function in the discrete 
approximation 

(e—1)(1—m,) 
e9 _{(1 —r)/(1—re!®)} un’ 





$, (8) = 


where r = m,n(1+m,n)-}, and this is the characteristic function of Z,, 
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where Z, is measured in units of n-!. To obtain a limiting distribution 
we therefore define a new variate Z; = Z,n-! whose characteristic 
function is y (4 (em _1)(1—m,) 

Dp - = ; - . 
"\n ein _f(] —r) (1—re®9in)\tin 


Letting n tend to infinity, where r = m,n(1+m,n)-}, we obtain 


(6) im 4, (2) = ; (l—m,)i6 ies (8) 
nao  \M 10 +-log(1—m, #6) 





and this is the characteristic function of the limiting distribution of 
Z, = Z,n~. In particular we find the moments of Z; to be 








i m3 

ws 2(1—m,)’ 

a. ee ae 

y (3(1—m,) 4(1—m,)?} 

2 el! m4 oe m> a me) 

Ms "40 —m,) * 30—m,? * 81 —m,)3)’ 

- oa { m3 m$ mj : me) 





\5(1—m,) 90 —m,)* a —m,)>~ 8(1 —m,)*} 
We can rewrite (8) as 
(1—m, )log(1—m, #@) 


(@) + 1 pe —_ — 1—m,+ ®, (8). 
) ™ 10+ log(1—m, 7@) . (0) 





Now a distribution function has a probability density if its characteristic 
function ¢(@) is such that 


x 
| 4(0)e~*° dé 
x 

exists as a principal value, and, when x + 0, this is clearly true for ®,(4), 
which is m, times the characteristic function of a continuous distribution. 
Thus 1—m, is the concentration of probability at the origin, as was 
suspected from the discrete case. F(Z), the cumulative distribution 
function of Z,, has a jump of magnitude 1—m, at the origin and is 
otherwise absolutely continuous. F(Z) for Z > 0 is therefore as given 
by Doob [(1) 39] 


F(Z)—3}{F(0+0)— F(0—0)} = F(Z)—4(1—m,) 


1 (1—m,)(1—e-*®2) 
27 } {i0+log(1—m, i8)} 





« 


which I have been unable to evaluate. 
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The above theory assumes that m, < 1 and considers the probability 
distribution at the bottom of a dam of indefinitely great capacity. Putting 
M = | in (1) and starting from the top of the dam we could set up an 
infinite set of equations for the distribution near the top of the dam 
which will have a meaningful solution when m, > 1, but the matrix of 
coefficients is then nearly upper-triangular instead of nearly lower- 
triangular as before, and the generating function technique does not 
work. 

Consider, however, the case of the finite dam. If we put M = 1 in 
(1), we see that P, is uniquely defined in terms of P,, P, in terms of P, 
and F, and so on. Let the solution of these equations be P4,..., Pe _,. 
Then P;/P, = P;/P, for the infinite dam. It follows that, if {P,, P,,...} 
is the solution of (4), the solution of (1) is given by 
, P, a” ‘ 

P; = B+P+..+h, (¢ = 0,..., K—1). (9) 
(Notice that (5) no longer holds.) This means that we can solve the 
discrete equations for the finite dam whenever m, < 1. But suppose 





_* 


that m, > 1 and consider the expansion of 


2-1 
z—p(z) 


= Apt 4,2+... 


in a power series. This always converges for sufficiently small z if py > 0, 
and the coefficients a; will satisfy the first K equations of (1), with M = 1, 
when substituted for the P;. The probabilities P; are then given by 


> a; 
ne Ag+... +dg_, 
whether or not m, < 1. 
By proceeding to the limit we see that in the case m, < 1 we can 
solve the continuous problem for a finite dam. The concentration of 
probability at the origin is 


1] —e-t9K -1 
jt |. dé}, 
J 10+log(1—m, i@) 





and the new cumulative distribution function F,(Z) is given by 
F(Z) = {F,(0+0)—F,(0—0)} 


2) 0) 


, | —e-i0z aol ( me [ | —e-i0K 
—_— +7 - po 1+7 2 5 dé bd 
\" ” i0+log(1—m,i6) | ° Bas . J 10+log(1—m, 76) J 








— 2 
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It seems probable that this formula remains valid when m, > 1, but I 
have not been able to prove this. 
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THE CHARACTERISTIC ROOTS OF THE 
PRODUCT OF MATRICES 


By N. A. KHAN (Aligarh) 
[Received 13 February 1956] 


1. Introduction 


DvuRING the last fifty or sixty years several authorst have given limits 
to the characteristic roots of an arbitrary matrix A = (a,;). In short, 
Browne [(3) 263] showed in 1930 that 


c(A)| < (R+T), (1) 
where 
R= max R&; = max ¥ 4a;,,, T = maxT; = max 9 a;;\, 
: . _ J . 5 ‘ _ t 
(%) (i) j ) Y) t 
and c(A) denotes a characteristic root of A. This result of Browne’s was 


improved by Parker (5) in 1937 who showed that 


e(A)| < max{}(R;+7;)}. (2) 
(a) 
The inequality (1) was further sharpened by Farnell (4) in 1944 when 


he obtained 
e(A)| < (RT)}, (3) 


and was still further improved by Barankin (1) in 1945, who proved 
that 


c(A)| < max(R;7;)!. (4) 
It is the purpose of this paper to give an alternative proof of (3) and 
to extend it to the characteristic roots of the product of two n-square 
matrices A = (a,;) and B= (b,;) whose elements lie in the field of 
complex numbers. In what follows the sum of the absolute values of 
the elements in the ith row of A will be denoted by R,(A), the sum of 
the absolute values of the elements of A in the jth column by 7;(A), 
and R(A) and 7(A) will stand respectively for the maxima of the R,(A) 
and of the 7;(A). 


+ References are given in (3) and (6). 


Quart. J. Math. Oxford (2), 7 (1956), 138-43. 
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2. The theorem given by Farnell may be stated as follows: 
FARNELL’S THEOREM. Let R, = '* a,,|, 7 = > |a,,|, 2 = max(R,), 

and T = max(T.). Then \c(A)! < (RT). ; 


To prove this let M<NM<... <2 
be the (necessarily real) characteristic roots of the Hermitian matrix AA’. 


Then there exists a unitary matrix U = (u,;), such that 


AA’ = UD(A2)U", 


where D(A*) = diag(A?, A3,..., A2), 
i.e. D(*) = U'(AA')U, 
so that i? > Ui; Xpg Ugis 

r,3 


where «,, stands for the element of AA’ in the (r,s)th position. This 
equation may be put as 


where 2,, %,..., 2, is a set of complex numbers such that ¥ 2,%, = 1. 
r 


Let |x, &., so that > & = | 
Taking the absolute value of (5) and remembering that A,’s are all real, 
we have 2 <F |a,|é,¢, 
r8 
> x, (2+ 3) 
r8 


°M 
WS 

“MM 
t 

*™M 


AA’), 


r 


=> GRA 
r 


since, AA’ being Hermitian, 7(AA’) = R,(AA’). Suppose that R,(AA’) 


attains its maximum value for r = k; then 


< R,(A4’). 


d? 


t 
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But sid 2 - ots 
R,(AA ) si > es hes > Ay, a, Peet | > Axs A ns| 

8s 3s 3 ' 


S > Ags | |4y3| 7 > |Aus||42q|+---+ > lxs| lng! 
8 8 8 


= Ides > Ay, | + |azo| > (ya | +... + [Opn | > ‘din| (8) 
= |ay,|T,(A)+ \aye|T,(A)+...+ |a,, |T),(A) 
< T(A)R,(A) 
< R(A)T(A), 
so that A? < R(A)T(A) (i = 1, 2,..., m). 
But, by a well-known result of Browne's [(3) 263] 
c(A)é(A) < ¢max(AA’) = dj, 
i.e. c(A)|? < A2 < R(A)T(A). 
Accordingly c(A)| < {R(A)T(A)}}. 


I now prove the following theorem which gives the upper bound for 
the modulus of an arbitrary characteristic root of the product of two 
n-square matrices A and B. 


THEOREM 1. Let c(AB) denote an arbitrary characteristic root of AB. 


Then c(AB)| < {R(A)T(A)R(B)T(B)}". 
Proof. Let #<AB<.. < 


A2 be the (necessarily real) characteristic 
roots of the Hermitian metrix 4A’. Then there exists a unitary matrix 
where D(A?) = diag(Aj, A3,..., A2). 
Similarly there exists a unitary matrix Q@ = (q,;) such that 
D(u?) = Q’( BB’)Q, (8) 
where D(p?) = diag(y3, u3,..., 42) 
and pw? < pz <... < p2 are the (necessarily real) characteristic roots 
of BB’. 
From (7) and (8) we obtain 
D(%u2) — PAA’ PQ’ BB'Q, 
i.e. AF pi ae > Pri Xs Pst Lut Pur Ieir 


7,8,t,U,v 
where «,, and §,,,, denote the (r,s)th and (u,v)th elements of AA’ and 


BB’ respectively. 
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Since P and Q are Hermitian matrices, the above equation may be 
written as 2 2 = . 
dj a Soe > pst, XB Is Yr (9) 
r,8,t 
where 2, %»,..., Z, aNd Y;, Yo,-.., Y, are sets of complex numbers, not 
all zero, such that ¥ z,%, = > y,9, = 1. 
Ss 


~ 


Let |x, é.and y, n,, 80 that 
F@=F 2 = 1. 
Tr s 


On taking the absolute values of (9), we have 


AR ME < Dd laps |€-F5\Bar! Ms 1 


Te 


wtA 
= 
‘A> 


: Ss > Bu Ns 
8,f 


i< 2 2 2 2)\ 
} j > Xrs (f; T &3) » 3 Bu (5 wm Ur )j 
r,s 8, 


4{ SRA) + SANAT)| |S of R(BB) +5 of NB) 
> RB, 1A’)€2) > R,( BB’)x}}, 





since R,(AA’) = T,(AA’) and R,( BB’) = 7, BB’), 
R,(AA)R( BB’), 
if R(AA’) and R, BB’) attain their maximum values for r = k and 
s = | respectively. 
But, as proved in (6), 
R,(AA’) < R(A)T(A). 
Similarly R,{ BB’) < R(B)T(B). 
Thus we have 
A? yu? < R(A)T(A)R(B)T(B), which is true fori = 1, 2...., n. 
By a recently proved result given by Roy (7), 
c(A B)e(AB) < Cmax(AA’ )emax(BB’) = 22 p?2. 
Hence, c(A B)? < R(A)T(A)R(B)T(B). (10) 
This completes the proof of the theorem. It may be observed that 
by putting B = J in (10) we get the inequality (3). 


I now establish another theorem which, in general, gives a better 
upper bound for |c(A B)} than that given by (10). 
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THEOREM 2. Let c(AB) denote an arbitrary characteristic root of AB. 


Then c(AB)| < min{ R(A)R(B), T(A)T(B)}, 


= 


i.e. the absolute value of an arbitrary characteristic root of AB does not 
exceed the smaller of the numbers R(A)R(B) and T(A)T(B). 


Proof. If A denotes an arbitrary characteristic root of 
AB = | ¥a;,6,,;), 
r 


there exists a set of complex numbers 2), £9,..., 2,,, not all zero, such that 


dx, = La, (¢ = 1, 2,..., n) (11) 


7 


and a set of complex numbers y,, Ys,..., y,, not all zero, such that 


Ay; = rp a. b (2 = :. | n). 


jr ri. 
7 


Let x; = &; and let &, be the greatest of the €;; then from 
Ay, = >> S Mey bp j 


irl Apje vj 


it follows that Ax,| = | > a, 6,;2;|, 
‘jr 


7 t<| 
1.€. A ox <> Uy b isi i> gr b,j (Ek: 


_ r).> 
3° 3.7 


But, since €, + 0, we have 
' 


AV < > Dp, by ai .) Ap, brs 7 sa > ky Brn 
r r 


r 


= (Ay > bys T |Ez2 p bo. Tee tT |Bpn > bus 
8s 8 8 


= 'd,,|R,(B)+ \a,. R,(B)+...+ |a,,,|R,(B) 
< R,(A)R(B), 
and hence lA! < R(A)R(B). (13) 
Now let \y;| = 7; and let », be the greatest of the »;, so that 
Ay, = > Arby y;- 
ir 
Taking the absolute value of the above equation and proceeding as 
we did in establishing (13), we obtain 
A! < T(A)T(B). (14) 


Combining (13) and (14) we establish the theorem. 
In particular, if we put B = J in the result of Theorem 2, we have 


c(A)| < min{ R(A), T(A)}, 
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a result due to Frobenius [(2) 387] for a matrix A with positive elements, 
und due to Brauer (2) for an arbitrary matrix A. 

It may be remarked that the limits given in Theorem 2 are, in general, 
better than those given by Theorem 1. However, when A and B are such 
= R(A)R(B) = T(A)T(B), 
then the two theorems give the same upper bound for |c(A B)). In par- 
ticular, if A and B are Hermitian matrices, then R(A) = 7(A) and 
R(B) T(B), so that c(AB)| < R(A)R(B). Thus we have 


THEOREM 3. Let A and B be two n-square Hermitian matrices and c(A B) 
stand for an arbitrary characteristic root of AB. Then 


c(AB)| < R(A)R(B). 


REFERENCES 

1. E. W. Barankin, ‘Bounds for the characteristic roots of a matrix’, Bull. 
American Math. Soc. 51 (1945), 767-70. 

2. A. Brauer, ‘Limits for the characteristic roots of a matrix’, Duke Math. J. 
13 (1946), 227-95. 

3. E. T. Browne, * Limits to the characteristic roots of a matrix’, American Math. 
Monthly, 46 (1930), 252-65. 

4. A. B. Farnell, ‘ Limits for the characteristic roots of a matrix’, Bull. American 
Math. Soc. 5041944), 789-94. 

5. W. V. Parker, ‘The characteristic roots of a matrix’, Duke Math. J. 3 (1937), 
484-7. 

6. ‘Characteristic roots and field of values of a matrix’, Bull. American 
Math. Soc. 57 (1951), 103-8. 

7. S. N. Roy, ‘A useful theorem in matrix theory’, Proc. American Math. Soc. 

5 (1954), 635-8. 








CENTRALLY BIASED DISCRETE RANDOM 
WALK 
By J. GILLIS (Rehovoth, Israel) 


[Received 28 February 1956] 


1. Introduction 
WE denote by m the general lattice point (m,, mg,..., mq) in a d-dimen- 
sional lattice and by E; the unit vector parallel to the positive direction 
of the ith axis. We now consider a random walk starting at the origin 
and such that the only steps permitted are of the type m > m+ E,, with 
respective probabilities P:(m), Q;(m) (¢ = 1, 2,...,d). A recurrent point is 
defined as one through which the walk will, with probability 1, pass an 
infinite number of times. The main purpose of this paper is to prove 
Theorem 3 below. However, the proof will require two preliminary 
results which it will be convenient to state separately as Theorems | and 2. 
THEOREM |. Jn a discrete random walk on a one-dimensional lattice let 
p;,; denote the probability of a step from lattice point i toj. Suppose further 


that Poa = Po-1 = 3: 


Pes = He) +1, +2....), 
Pisa = 2(1+e/t) 
Pi3 = 9 when \i—j| +1, 
where \e| <1. 
Let the walker start from the origin, let B, denote the probability that he 
will return to it, not necessarily for the first time, after r steps, and write 


By = 1. Then, for |z| < 1, 
¥ B.2t = F(je+1, Je+-4; 1s2*)/F(de, fe+4; 152), (1.5) 
r=0 


where F is the elementary hypergeometric function. 


THEOREM 2. Let {u,,} denote an infinite sequence of non-negative real 


numbers such that, as n tends to infinity, 
r n—1 
u (1.6) 


r 
r=0 


where X is constant. Then, given 6 > 0, we can find K = K(A, @) such that 


awit <u, < Kn)-1+8 (1.7) 


for all sufficiently large n. 


Quart. J. Math. Oxford (2), 7 (1956), 144-52. 
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THEOREM 3. Consider a random walk on a d-dimensional lattice which 
starts at the origin and is such that 

(i) ifm, 40, Pm) = (1—e/m,)/2d ee 

(¢ = 1}, 2...., d), (1.8) 
Q,(m) = (1+ .«/m,)/2d 
while 

(ii) ifm,;=0, P(m) = Q,(m) = (2d) (i= 1,2,...,d), (1.9) 
where «| <1. Then the origin will or will not be a recurrent point according 
” d—2 _d—2 

ee ry ~ or i se 

Remark. The case « = 4—d~ has not yet been elucidated in general. 
However, when d = 1 it is known that e — 4 makes the origin a recur- 
rent point [cf. (2) and (3), or the result can also be easily deduced from 
(1.5)]. Again, when d = 2, « = 0 corresponds to unbiased random walk 
in a plane, and it is well known (4) that in this case the origin is a recurrent 
point. It seems a plausible conjecture that, for all d, «e = }—d-' makes 
the origin a recurrent point. 

Under the conditions of Theorem 3 every point is accessible from every 
other point. It follows (4) that all points are recurrent if and only if the 
origin is. 

We note that, by Theorem 3, « = } will make the origin a recurrent 
point in a lattice of any finite number of dimensions. 


2. Proof of Theorem 1 

For all n > 1, let a(m,n) denote the probability of arriving at the 
lattice point m on the nth step. It is convenient to define a(m, 0) to be 1 
when m 0 and zero otherwise. Then, for x > 1 


> 


x(M,N) = Pm-1.m(m—1,n—1)-+p a(m+1,n—1). (2.1) 


m+1,m 


Multiply both sides of (2.1) by zz" and sum for —0%0 <m<ao,n >1. 
This gives 


A(z,z)— > a(m,0)a™ = xzP(x,z) + - Q(z, 2), (2.2) 
m=-—® zx 
where 
ao 2 
A(z,z)= YY } alm,n)x2", (2.3) 
m=—-x2 n=0 
P(z,z)= > > paws: alm, 22%", (2.4) 
m=—x n=0 
a @ 
Q(z,z)= SF > Dam-12(m, 2x2". (2.5) 


m=-—D 


0 
L 
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It follows from the definition of a(m,0) that (2.2) can be written 
rP+Q/x = (A—1)/z, (2.6) 
where, for simplicity, the arguments of the generating functions A, P, Q 
are omitted. Again we know that, for all i, p; ;_;+p;;,,; = 1 and hence 
P+Q= A. (2.7) 
From (2.6) and (2.7), 
P = {x(A—1)—zA}/2(2*—1). 
Now, by (1.1), (1.2), (2.4), 
P = 4A—Jfe SWE) gman, 


_ m 


where >’ is taken over all m except m = 0 and all n > 0. Hence 
oP 1@:z 
SF is an be S” a(m, n)x™—12" 
Cx 202 *™ 
1éA 
= =~ }ex-(A—B), 
= OL 
where B(x) = > a(0,n)z" = > B, 2’. (2.11) 
n=0 r=0 
But, by (2.8), 
2(x?7— 122P = A[2zx—(2z?+1)] = (a?—1)(~a—z)+(2?+1). (2.12) 
cx 


Cx 


Comparison of (2.10) and (2.12) leads to an ordinary differential equation 
for A as a function of x. To solve this equation we write 


x= e'P, (2.13) 
A(x,z) = E(¢,2), (2.14) 


and the equation becomes 


(l—z cos 4) + [z(1—e)sin d—(1—zcos¢)cot d]E 


c 


= —cot d—«zB(z)sin ¢. 


This is easily solved by means of an integrating factor, yielding 
E(¢,z) = (1—z cos ¢)-!+ ez sin (1 —z cos ¢)-!** i (1—z cos @)-!-« d6— 
‘ 0 
—ezB(z)sin ¢(1—z cos ¢)-1+€ | (1—zcos 6)-€ dé + 
0 
+sin d(1—zcos¢)-!**K(e, z), 


where K depends only on the parameters indicated. 
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It is clear from the symmetry of the conditions that 
x(m,n) = a(—m, n), 


and so E(¢,z) is an even function of ¢. The first three terms on the right 
of (2.16) are also even while the last is odd, and hence 


KE = 0. (2.17) 
Now E(¢,z) = > a(0,n)z"+2 > + af x(m, n)z” cos md, 
n=0 m=1 n=0 
and so B(z) xf E(¢, z) dd. (2.18) 


The integral in (2.18) can be evaluated from (2.16) and (2.17) if we 
integrate by parts the second and third terms on the right of (2.16). We 
then get 
27 
{ E(d,z) dd 
0 
an 


(1—z)£ | (1—zcos @)-!-€ d8— 
0 


o — B(z)(1—z)¢* ( (l1—zcos @)-€ dé +27B(z). (2.19) 
0 
From (2.19) and (2.18) we get 
B(z) [ sie an + ao) (1—z cos @)-€ dé, (2.20) 


0 
and this is precisely (1.5) as may be seen by expanding the integrands 
in (2.20) in powers of z and integrating term by term. 


3. Proof of Theorem 2 


Given 6 > 0 we can find ny = n,(@) such that, for all n > np, 
A+86 n—1 
u, << — U,. (3.1) 
? 
r=0 


} by 


nj 


We now define a new sequence {v 
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It is clear from (3.1), (3.2), (3.3) that 
v. >u, forall r. (3.4) 


r 4 
We can prove inductively, using (3.3), that v!’” is bounded and so the 
series > v,2’ has a positive radius of convergence. Within the circle of 
convergence we write 
(3.5) 
It follows from (3.3) and (3.5) that, for sufficiently small ‘z , 
(1—z)-*V(z) = v9 + (Up +0,)2+ (U9 +0, +0,2)2?-+4 
(A+6)-? $ rv, 2’-1+-A,(z), 
r=] 

where A,(z) is a polynomial in z of degree not exceeding mp, i.e. 

dV 2X+6,, 

a aeneenes PEE) as eee. (3.7) 

dz l1-z 
where A, is also a polynomial of degree not exceeding n,. Integrating 


this equation we get 


V(z) = vo(1—z) 4" + (1—z)-*-” [ (l—u)+"(ag+a,u+...+a,,u™) du, 


0 
(3.8) 
where the a;’s are constants. Now it is easily seen by successive integra- 
tion by parts that 


(1—z)-4-? i u'(1—uy\+? du = RAz)+S(1—z)?~, (3.9) 
0 
where R; is a polynomial of degree i+1 and S; is a constant. Hence, 
by (3.8), V(z) = K,(1—z)?+ Lie), (3.10) 
where K;, is a constant and L a polynomial of degree not exceeding ny+ 1. 
It follows that, for n > n,+1, 
- (—A—84\ K, n4+9-1 ; 
v, = K Pp eglieremmae « 3.11 
( n (A+ 9) \ 
From (3.4) and (3.11) we immediately deduce one of the inequalities 
(1.7), and the proof of the other inequality is completely analogous. 


4. Proof of Theorem 3 
Let us denote by 8 the probability that the walker wil] return to 
the origin, not necessarily for the first time, after r steps. Then it is 
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known [ef. (1) 244] that a necessary and sufficient condition for the 
origin to be a recurrent point is that 

> Aa 

) 

> Be 

r=0 
diverges. We define BS” to be 1 and we recall that, in the notation of 
Theorem 1, 

BY = B, = a(O,r). 

The numbers 8 depend, of course, on ¢. In the interest of clarity I 
shall formulate the proof as a series of lemmas. 


LEMMA 1. 


B(z) = > B27 = 1+ 5 (log F(4e, 4¢+ 4; 1; 2°)]. (4.1) 
« dz 


r=0 


Proof. Writing u 2? we have 


F(a+1,6;c:u)—F(a,b;c;u) = . F(a+1,6+1;¢e+1;u) 
- 


oem. qd F(a,b;c;u). (4.2) 
hbdu 
Hence, by (1.5), 
B(z) ] 2zt _d flog F(4e, 4 1; 1;2%)] 
(2 ae — — y S¢€, S€-- 8,452 
— 


zd - 
: I+- 7,108 Fide, de+45 1524). 


€ de 
LemMa 2. B(z) satisfies the Riccati equation 


dB 2 + 2e)z Zz 
* — Spr [=+O°2*|p_s_ ;: (4.3) 





Proof. This follows immediately on substitution from (4.1) into the 
differential equation for F. 
LemMaA 3. As N tends to infinity, 
g e+ > Q 
Pen “™~ N ra Pp 
=0 


(4.4) 





2i° 


It is clear from the nature of the problem that £, is zero for odd r. 
We write, then, 
x 
3(z) = S y,2** (where y; = B3;), (4.5) 
i= 
insert this series in (4.3), and compare coefficients of the various powers 
of z. This gives 


N N-1 
2Nyx = —€ 2 viyy-at ey t(1+ 2) ¥ ni. (4.6) 
i= 
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Suppose that [Ze] < 1. (4.7) 


We see by letting z tend to 1 in (2.20) that 


= « 
> yy = +0. (4.8) 
N=0 
We next wish to show that 
lim yy = 0. (4.9) 
N--« 


Consider first the case « < 0. Then 


N-1 N-1 
2Nyyx (1+2e) ¥ y,—e 3 vivy-i—e, by (4.6), since y, = 1, 
t=0 1=0 
N-1 
(l+e) ¥ y;—e, since 0 < y; < 1 for all i, 
t=0 
N-1 
< (1+4e) ¥ y;, for all sufficiently large V, by (4.8), 
i=0 
e+2 : =! 
i.e. Yn <- es S Yi (4.10) 
t=0 


and so, by an argument precisely similar to that used in the proof of 
Theorem 2, , 
YN O(N-!-«4) = o(1). (4.11) 


When « > 0, we have, from (4.6), 


2Nyy < (1+2e) ¥ y;, (4.12) 
i=0 
and so, as before, 
Vy O(N*-!) = o(1)_ by (4.7). (4.13) 


Now take any arbitrary fixed positive 5 and choose m, so that y; < 8 
for alli > m,. Then, for V > mp, 


N mo N 
Yvivv-i<Sut+ > vivy-w simee y; < 1 for all, 
i=0 i=0 i=me+1 - 
mo N 
7 th hs & 
ho Ft 0 -_ YN-i 
i=0 i=me+1 
mo N 
b 3 vito > Vie (4.14) 
i=0 i=0 
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From (4.6), (4.8), (4.16) we deduce that, for all sufficiently large NV, 





N-1 



























¥ 2% (4.17) 


which is (4.4). 
LEMMA 4. Given any positive 0 we can find K, = K,(@) such that, for 
all N, K-1Net0 < y. < K,Ne4, (4.18) 
Proof. This is an immediate consequence of Lemma 3 and Theorem 2. 


LemMA 5. Given X > 0,A > 0, p > 0, we can find K,, K, depending 
only on A, p such that 
n 
0< K, < n“(X+1p5 (3) wipe < Ky. (4.19) 
imo \*" 

Proof. This follows by considering the largest term of the sum in 
(4.19) and estimating this and neighbouring terms by Stirling’s formula. 
The proof is closely similar to that of the normal approximation to the 
Bernoulli distribution and the details will be omitted here. 

LEMMA 6. Given any positive yp and a positive integer d, we can find 
K,, K,, depending only on p, d, such that 

d 
0< K,< dr - (2r)! [] [(2n,)!n4]-? < Ky. = (4.20) 
ny +...+Ng=T i=1 

The sum in (4.20) is to be taken over all sets of positive integers 
Ny, Ny,..., Ng Which satisfy the indicated condition. We prove the result 
by induction over d. The lemma is clearly true when d = 1. In the 
general case write ? 

wb(d,r, w) d-2rydu > YT L [(2 n;) ! mH]- (4.21) 
Nit. TNa= = 


and assume the lemma for d—1. Then 


r—1 
ub(d, 1, w) d-2rydp > [(2n 7)! n4]-(d—1)"-"4(r— Nq) 4-4 (2r)! x 


na=1 
x [(2r—2n,)!]b(d—1, r—ng, 2) 
r—i 


oo, 
< K.(d—1, p)d-**rt at i—4-De(r —i)-#(d—1)*, 
I 2 


by the inductive hypothesis, 
< K,(d—1,)K,[u,(d—1)u], by Lemma 5, 
- K,(d,u) (say). (4.22 


The inductive proof of the lower bound for y is similar. 
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We can now complete the proof of Theorem 3. Under the conditions 
of the theorem let 8 denote the probability that the rth step will end 
at the origin. As before, 8 is zero for odd r, and we write yO = p. 
Then it follows from elementary probability considerations that 


: » V)! 
y@ = d-2 9 a vy ¥,. (4.28 
. 4 : (2n,)! (2mg)!...(2mq)t 27 Yn | 





ot 
Ny +N... + ng=N 
Suppose that « < (d—2)/2d and choose @ such that 
e+0 < (d—2)/2d. (4.24) 


rr . . . a . . a s . €—3+0 
rhen, for each i, Yn, < Ay nf : 


by Lemma 4. Hence, by Lemma 6, 

2) — Kad g-2N = IN)! a In .\!]-1 4 =—e—<¢—} 

ve < Kid a (2N)! TT (22) TT 2; 

y+... ng=N i=1 i=1 
< K{K,(2N)-#4-<-, 
and so, by (4.24), ¥ y@ converges. Similarly we can prove that > y\? 
diverges if e > (d—2)/2d. This proves the theorem. 
It remains to consider the restriction imposed by (4.7). The condition 

e < } clearly changes nothing since we have shown that, for all finite d, 
the critical value of ¢ is in fact less than }. For a similar reason the 
condition « > —} can affect only the instance d = 1, where Theorem 3 
is in any case known to hold (cf. the remarks at the end of § 1). 
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A SQUARE AND A PRODUCT OF 
HYPERGEOMETRIC FUNCTIONS 
By H. F. SANDHAM (Dublin) 

[Received 16 March 1956] 


| sHow that Clausen’s (1) and Orr’s (3) theorems 


\3 F( 2a, 28, a+ 8; ). (1.1) 


1. 


Zz 
x+B +4, 2a+ 28 


xB; _\ pf/t—a,4—B; \ oft. $+-0—B, $—a+8B; || » 
pea) F gap) =F crpth deep 7) 


-although 70 years separates their diseovery—can be deduced directly 
from one another on our recalling the classical theorem, which for hyper- 


geometric functions of order 3 is that 


r(" a ty *) 
aa p(4 + 1—d, b+-1—d, c+1—d; 
P| 2—d,e+1—d i 


\ 


rie p(tt+l—e, b+1—e,c+1l—e; a 
, d+-l1—e, 2—e 


satisfy the same linear differential equation of the third order with 
polynomial coefficients. 
Also I deduce Orr’s theorem from the Dougall-Ramanujan formula (2) 


Cc, d, é, —m; 1) 





F(* 1+ 3a, b, 
la, 1+a—b, 1+a—c, 1+a—d, 1+a—e, 1+a+m 
(1+a),,(1+a—b—c),,(1+a—b—d),,(1+a—c—d),, 
(1+a—b),,(1+a—c),,(1+a—d),,(1+a—b—c—d),,’ 
where 1+2a = 6+c+d+e—m: 


and thus show that these two most interesting theorems on products 
of hypergeometric functions can be deduced from well-known results, 


without calculation. 
Quart. J. Math. Oxford (2), 7 (1956), 153-4. 
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2. On expanding the hypergeometric functions in powers of x and 
equating coefficients of x" we see that Orr’s theorem is equivalent to 
rk eee x; B, yn) = Walk re Path ath 
B a—n, a+}—n, xr+B+4 (4—2)n($—B)n(a+B+4),’ 


which is a particular case of the Dougall-Ramanujan formula. 





3. We next note that 


Q. \ : = 1 2p. 
A= P| », B; z), B=2 =-OR(*  2—Ps 
x+-B+4 a 
Rp > > x Pp } 


satisfy the same linear differential equation of the second order with 
polynomial coefficients and that therefore the lowest-order linear differ- 
ential equation with polynomial coefficients satisfied by their product 


has only solutions 42 AB. RB 


and is therefore of the third order. Orr’s theorem tells us that one 


solution of this is 


on 28 p(/ 1 —2a, 1—28, l1—a—B: 
c) x} se BF é *). 


$—a—f, 2—2a—28 


92° 


These are easily identified with A? and B*, though the latter identifi- 
cation merely yields Clausen’s theorem once again, with a, 8 replaced 
by 3—a, ? 
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CHARACTERIZING PROPERTIES OF 
STATISTICAL DISTRIBUTIONS 


By J. G. MAULDON (Ozford) 
Received 21 March 1956] 


Ir x and y are independent variates and x+y is known to have a normal 
distribution, it is well known [Cramér (1)] that z and y must each be 
distributed normally. It follows that, if the distribution of the mean 
of a sample from a population is normal, then the population distribution 
must also be normal. This suggests the general problem of finding out 
when the distribution of a sampling statistic determines that of the 
population. 
PROBLEM 1. Suppose that 2,, %,..., 2, are independently distributed 
in the same distribution defined by 
dF = f(x;)dz; (—2o <2,< 0; §= Il, 2...., n). (1) 
Let %, s > 0, t be defined by the equations 
nz ¥ x.. (n—1)s? > (x;,—Z)?, t ENn’s, (2) 
and suppose that ¢ is distributed in Student’s ¢-distribution (7). Can we 
then assert that the distribution defined by f(x) is the normal one? This 
will be shown to be false for n = 2 and is probably false for higher values 
of n. It has been conjectured by P. A. P. Moran and the author that the 
conclusion is true if the hypothesis is assumed true for all n, but this we 
have not been able to prove. 
T he Case i 2. 
#=H}a,t2,), 8=2-4\z,—2,), t= v2. (3 
| 1 2 1 2 
The joint distributions of (71,22), (%, 8), (t,8) are 
dF = f(z,)f(z.)dz,dz, (—© < 2, 2%, < ®), (4) 
dF = 2!f(g+2-'s)f(g—2-1s)dids (—-w << F< D;0<8s< @), (5) 


dF 2f{(t+-1)s/v2}f{(t—1)s/s2}sdsdt (—w2w <t<0;0<8< o). 


(6) 
If the marginal distribution of t is 
dt 
dF > > —-D <= t <= ©), 7 
mise) | (7) 


Quart. J. Math. Oxford (2), 7 (1956), 155-60. 
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 * 
then 2r | f{(t+1)s/v2}f{(t—1)s/v2}s ds = 


« 


0 


l 
1+-# 


for all values of ?. 
We now write 
u 4s*(1-+1)?, a = (t—1)?/(t+1)?, 

g,(u) vaf(\ut}), J(u) = Vaf(—|ut)). 
Then, supposing first that t > 1, we find that (8) becomes 

| 9,(u)g,(au) du = (l+a)"?* (O<a< 1). (11) 

0 
Equation (11) still holds if a > 1, as we see by writing au = v and 
a~! = b, and a continuity argument for the case a = 1 gives (12) below. 
By considering the cases —1 <t< 1, t << —1, we find the further 
results (13), (14). 


( g,(u)g,(au) du = (1+a)-! (a > 0), (12) 
i 
| g,(u)g.(au) du = (1l+a)-! (a > 0), (13) 
0 


( g.(u)g(au) du = (l+a)-! (a> 0). (14) 
0 

Adding equations (12) and (14) and subtracting twice (13) gives, 

when a= l, es 
[ {9,(u)—g.(u)}* du = 0, (15) 

0 
from which it follows that g,(u) = g.(u) almost everywhere and there- 
fore that f(x) is effectively an even function (i.e. ignoring a possible set 
of measure zero). Hence the most general solution of equation (8) must 


be of the form f(x) = a-4g(x2) (-0w <4 <0), (16) 


where [ g(u)g(au) du = (1l+a)-! (a> 0). (17) 
0 
This integral equation (17) has been studied by Goodspeed (2) and 
has an infinity of solutions of which the following three have been 
taken from his paper. 


(a) g(u)=e-*, f(z) =a (—wo <2 < oo). 
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This gives the normal or Gaussian distribution, and the result holds for 


any change of scale. 


}il+u* 


‘ 9 
(b) g(u) (=| ; f(x) _—_—— (—70o < 4% < O&). 
A \n TW Lx 
This is the simplest non-normal distribution with the required property. 
q(u) f(x) a (—o < xz < oO). 
t a(i-+-2Z 
This gives a bi-modal distribution which has the improper mean zero but, 


like the Cauchy distribution, has no proper mean. 


PropLeM 2. If x and y are independent variates distributed in 


identical y?-distributions, so that 
dF = {T(a)}-*(ay)*"e-** dady (0<2,y<w;a>0), (18) 


then their ratio v = y/x has the particular F-distribution (20), but this 
property does not characterize the y?-distribution uniquely. 
Let the joint distribution of x and y be 


dF =f(x)f(y)dady (O<2,y< @). 
Let v = y/x, and let the distribution of v be 


dF - nay” (20) 


Then (19) leads to an alternative form for the distribution of v, viz. 


dF = dv f(a) f(vx)x dx (0<v< oo), 


0 


and this is identical with (20) if and only if 


T'(2a)v4-" 
{ T'(a)}?( es y)24 





f(x)f(vx)x dx = (v > 0). 


. 


0 
Write [(a) f(x) = a4—19(x) 
and substitute this expression in (22). We find 


. 9 
g(x)g(va)x24-1 dx = Tae (v > 0) 
1 9;)2a 
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and Goodspeed [(2), Theorem 5] has studied this equation also. The 
general solution is 


g(x) = x [ k,(2*)e-* dt, (25) 
0 


where k,(u)/u is any L? Watson kernel, so that [Titchmarsh (3), (8.5.3), 
(8.1.9)] 


T 
ko (2) = 5 him. a! Bal dt, (26) 


-—T 


where K($+2t)K(4—it) = 1. 


I shall confine attention to functions &(s) satisfying the following 
three further conditions 


RK(4-+it)| = o(1) (t real), 
and S(4+it) satisfies in the neighbourhood of t = 0 conditions for 
the convergence of a Fourier series. (30) 
Then—the proof is given after (34)—we find that 


[ x-tk,(x) dx = 2 
a 


and hence, from (23), (25), (26), 


{ f(a) da = 1, (32) 
0 


and so the function f(x) gives a genuine probability distribution provided 
that f(x) > 0 for all z > 0. 

This is certainly true if k,(x) > 0, for example if 

l ] 
Rs)=1, klzy= ft &> 2) (33) 
\0 (0<2< 1), 

which leads to the y?-distribution (18). 

Other examples, which prove the statement at the beginning of this 
problem, are provided by (43) and by 


K(s) = | (E)reminazey k(x) = [E)a—cosa >0. (34) 
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Proof + of (31). 
























x x T 
rth, (x) dx = = l.i.m. | dx x-* [ RE +H) yaa dt 
} 20 T+0 , J 4-4 
x x71 -T 
x 
ROR [a R(3 + it) [« 1-8 de 
27 J 4—it . 
<a x-1 


dt 





a: vad xX-#_ xX 
ae 4—it —it 


é.. 
a7 


* 
—@ 





J t 


_ - Lf ‘SG it ) 2sin(¢log X) 5, 
2a ,— ; 


Using (30) and letting X — 0, we find that 


az-tk,(x) dx = sat = 2, by (28), 


. = 


0 


which completes the proof. 
Similarly 





: - - 9Q(1 9 
Hh, (x) dx ] iim 2K(4-+-it)sin{(t+ta)log X} dt. 
J 277 Xm | (4—1t)(t+ta) 
0 — 2 


and hence (by contour integration), provided that |a| < } and K(s) is 

regular for 4 x} <res < $+ al, 
és 

2% (x) dx = —=-—. (35) 


0 
PrRoBLEM 3. If x and y are independent variates and if v = y/x has 
the general F-distribution (37), then (as is intuitively plausible) it is still 
true that x and y need not be x?-variates. 
Suppose that 
dF = f,(x)f(y)dady (0 <2, y < @) (36) 


and that the marginal distribution of v = y/x is defined by 


7” net) dr (0<v<o@;a,b>0) (37) 
—— yee - , , in : 
[(a)P(b) (1+) “9.4, : 








+ This proof is due to Professor Titchmarsh. 
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Proceeding almost as before, let 
P(a)f,(x) = Ax*—y(x), (6) fa(y) = A*y’g(y). (38) 


Then we find 


x 


(a+b) 


[ ga)g(vx)e2t dx = 
0 


and this is satisfied by 


g(x) = 2 ( k, (t**°)e-= dé, 
0 
where k,(x) is given by (26). 
Using (35) with 2a +(a—b)/ (a+b) we get 


. } \ = 
f(x) dx AR{—— , foly) dy = AR =~ ) (41) 
a+b i a+b, 


0 
By virtue of (27), which can be written &(s)R(1—s) 1, we see that 
both these integrals take the value 1, provided that 


_@ \_ {of 4 \\" 49 
A a(— ; a 8(=)| j — 


. 


Hence f,(2) and f,(y) yield genuine probability distributions if, for 
example, k,(2) is as defined in (34). Another example is provided by 


l~+-1]} 
oe (43) 


= 


] 
R(s) = cot(4zs), b,(z) = —log | 
7 |z— 1} 
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